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ABSTRACT 

Complementary binary sequences were invented by Golay in an in- 
vestigation of infra-red multi-slit spectrometry. This dissertation 
formalizes the basic results obtained by Golay and develops new concepts 
and techniques for examining the characteristics of these special binary 
codes. This work has developed new understanding of the structure and 
methods for the decomposition of complementary sequences, 

Complementary sequences have the property of an infinite correla- 
tion peak to ambiguity ratio when detected with a matched filter, These 
binary sequences should find much application as pseudo-random noise 
modulation signals for both radar and communications systems. 

A discussion of the need for such sequences is included in the intro- 
duction and is followed by a state of the art description. An operations 
group on the sequences is formulated and the proofs of several theorems 
concerning the operations group are given in a rigorous manner, One 
reason for developing the operations group is the application to elimina- 
tion of redundance in the computer search for new codes. 

Several invariant properties of complementary codes are proved 
through the use of the Hamming distance concept, Many more invariant 
properties of the sequences are demonstrated through the introduction of 
a Hamming vector. The concept of a Hamming vector is extremely use- 
ful as a complementary code decomposition tool, A large number of 
theorems are proved to enhance its use in this field. 

Several computer searches for complementary sequences are 
described and the actual computer programs for the CDC 1604 are in- 


cluded in the Appendix, 
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CHAP? Ewa 
INTRODUCTION 


Search radars of today are required to obtain smaller and smaller 
targets at ever increasing ranges. This has led to many different tech- 
niques of summing radio frequency energy to obtain these long ranges, 
Complementary sequences are well adapted for coding the RF energy 
pulses to give a summing of the returning RF energy inthe echo. The 
present work is more concerned with the characteristics of complemen- 
tary codes suchas their operation groups, Hamming distances and 
Hamming vectors, than it is with their application as modulating signals 
for RF carriers, However, the introduction will develop the radar 
problem as a vehicle to emphasize the use of these codes in furore appli- 
cations in radar. Their application will probably be equally important 
in the communications field, The initial work in complementary sequen- 
ces was directed toward an infra-red application, but it is felt that the 
eventual application of these sequences will spread to many different 
fede. a 

The usual approach to the problem of increasing radar range is to 
look for ways of varying some of the parameters in the radar range 
equation which will lead to a greater maximum range. 

The radar range equation is: 

‘[P, GA_ SD 


R act eed 


me) | Peele 
min 


1/4 


where P. = peak transmitted power in watts, 
P_.. = minimum peak detectable signal in watts, 


min 


S = scattering cross section of target in units 
consistent with range, 


G = gain of transmitting antenna, 


A_ = effective area of receiving antenna in 
consistent units, 





D = acomposite loss factor for transmission lines, 
atmospheric losses, etc, 


The parameters under the control of the design engineer are AG, 
eer and P A. and G are fixed by various considerations such as 
antenna beamwidth, side lobe suppression, physical size of the antenna, 
etc. , and are not truly available as adjustable parameters beyond small 
variations, Therefore the two parameters for the designers to optimize 


are F . 
m 


and F.. 
in t 


ere is basically limited by three factors: Johnson noise in the 
input circuits, shot effect and other noises in the first tube, and cosmic 
noise picked up by the antenna, P ain approaches KBT_ in a perfect 
receiver, where K is Bolt zmann's constant, B is the bandwidth and a 
is the absolute space temperature, Fae is currently being reduced 
toward the above value through the utilization of parametric amplifiers, 
masers and other low noise devices, Another method of reducing Pee 
is through pulse integration, since noise voltage averaged over a period 
of time has a value approaching zero while the signal values are additive, 
However, in most cases this integration of pulses can also be considered 
as an increase intransmitted energy, It will be considered from that 
viewpoint in this paper, There are two basic types of integration, co- 
herent and noncoherent, Coherent integration gives a gain of approxi- 
mately N, where N is the number of pulses integrated and noncoherent 
integration gives a gain of approximately /N, Both of these values are 
compared to a single murse a The values are respectively the upper 
bound for coherent integration and the lower bound for noncoherent 
integration, 


There are basically three different integration time bases possible 


in a radar system. One is the antenna scan to antenna scan where the 





Operator or a computer, decides if the echoes from previous scans were 
signals or noise, The second method of integration 1s the interpulse 
method where the cathode ray tube, a delay line, or a computer sums up 
the hits within one antenna scan of the target. The third method is intra- 
pulse integration, or pulse compression, where different portions of 
one radar transmitted pulse are summed up, Complementary codes are 
used in one scheme of intrapulse integration, 

An examination of how the transmitted power or energy since 
BE, = P. A T, can be variedis now inorder, From this relation, it is 
seen that there are two basic methods for increasing the transmitted 
energy: either increase peak pulse power, or increase the pulse width. 
A common method of increasihg radiated energy has been to increase 
the peak power. This method has proven both expensive and wasteful as 
many radar transmitter tubes are peak power limited while only dissipa- 
ting a very small percentage of their average power capability, and 
radar modulators are taxed with a very similar problem. The alterna- 
tive did not look very promising at first, since to increase pulse duration 
normally reduces the high frequency content of a pulse, which in turn 
reduces range resolution capabilities, However, intrapulse modulation 
retains the high frequency content of a narrow pulse, while increasing 
the energy content to that of a pulse of long duration, Intrapulse modula- 
tion allows the summing, required in intrapulse integration, for 
detection in a radar receiver, 

There are two methods of intrapulse modulation, the digital method 
and the analog method, Although the emphasis in this paper is on 
digital modulation schemes - the analog system will now be briefly 
mentioned since it is a part of the overall picture, Figure 1.1 contains 
the graphs used in the explanation of the chirp radar elaine In 
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FIGURE 1,2 
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this analog method of intrapulse modulation the frequency of the outgoing 


pulse is swept linearly upward with time as is shownin Figure 1.1-b. 
This gives a large high frequency content to the pulse while at the same 
time maintaining a long pulse duration. Upon reception the signal is 
passed through a reversed time delay, delaying the low frequencies more 
than the highs, as is shown in Figure 1.l-d. This allows all the fre- 
quency components of energy to arrive at the output simultaneously, 
giving a large and narrow output pulse. The type of filter used in the 
delaying process is knownasa matched filter and is discussed in 
Appendix I. The procedure here is seen to be quite simple: first we 
generate energy components with various frequencies at different times, 
and then we delay all components so that they arrive at the output 
simultaneously. 

As a preliminary to the digital type intrapulse modulation consider 
a transmitted waveform suchasin Figure 1.2, This waveform con- 
sists of three pulses of unit length with spacing between them of one and 
three time units respectively. This waveform is then fed into a delay 
line in the receiver which has the delays shown in Figure 1,2-b, The 
output is either a one or a zero except at the exact match of signal to 
filter; the output is then three as is shown in Figure 1,2-c. The three 
pulses have therefore been summed as though they were one, at the 
time between 6 and 7 units. These pulses still contain all the high 
frequency components of the short pulse of unit length but have effec- 
tively three times the energy transmitted ina single pulse at the time 
of exact match. There are two features to be noted about this output, 
first it has been delayed in time and second there are secondary pulses 


at various undesired times which might be classified as coding noise, 
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Fig. 1,3 - Matched filter detection of a random code, 


Fig. 1.4 - Matched filter detection of an optimal noise code 





I will consider any undesired signal which might be mistaken for a true 
signal as noise. These two features will be discussed in some detail 
in the true digital intrapulse modulation schemes to be discussed later, 
The scheme just considered is somewhat inefficient since seven 
units of time are required to radiate three units of energy. If in the 
transmitted signal all the spaces with zero energy output in the pulse 
envelope had also had an output, seven units of energy would have been 
radiated in seven units of time. One might use some sort of frequency 
shift or phase modulation to accomplish this purpose, still utilizing of 
course the same output tube. Consider the same waveform as before 
with the units of energy coded as plus and minus by the modulator by 
means of phase shift modulation of either zero or 180°, The waveform 
is as shown in Figure 1.3-a. Upon reception this waveform is now fed 
into a delay line inthe receiver as shown in Figure 1,3-b. "A' repre- 
sents a straight-through amplifier, Figure 1.3-c shows the output of 
the delay line summer with respect to time. The output energy is now 
seven times that of a single pulse and has the frequency content, and 
therefore the range resolution, of one of the single pulses. A time delay 
in the output is again apparent, as are a number of noisy sub-peaks, 
This was a code picked at random and the results were good, There are 
certain classes of codes which are optimal, having even smaller noisy 
sub-peaks, These codes are of length 2-1 and have been studied 


5,6,7 They are knownas L codes or pseudo-random noise 


extensively. 
codes, An example of one of these pseudo-random noise codes is 
-+--+++, where the minuses and pluses indicate the phase of the RF 
signal, When this particular noise code is impressed upon its matched 
filter, Figure 1,4-b, the output of this summed delay line is as shown 


in Figure 1,4-c. It is to be noted that the maximum height of the 
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Matched filter detection of a complementary 
code of length 8, 


FIGURE 1,5 





"hash'' (sub-peaks) is one, with a main peak of seven, whereas in the 
previous example picked at random the hash had a height of two. 
Elaspas has made a study and design of a radar utilizing this form of 
modulation and matched filter detection, : He infers that codes of 
reasonable length should be used, For example, a code of length 127 
would be a reasonable length to obtain the benefits of a pseudo-random 
matched filter radar, An optimum code of this length has a maximum 
subpeak of 13, which gives a good coding signal-to-noise ratio, 

There is a way, however, to completely eliminate the hash, al- 
though this required the transmission of two separate sequences. Both 
of these sequences could be transmitted through one output tube with a 
moderate band width. These codes which eliminate the hash are known 
as complementary codes, 7 An example of a complementary code pair 
being detected by a matched filter pair is shown in Figure 1.5, The 
two codes being transmitted will be called A and B for convenience, 
with their respective matched filters also being labeled ina similar 
fashion. The receiver output, after the complementary codes are detec- 
ted by the matched filter and their outputs summed, is shown in Figure 
1.5-g. This demonstrates, as was mentioned before, that the hash level 
is zero, This lack of hash is indicative of an infinite coding signal-to- 
noise ratio. This particular feature is the prime advantage of complemen- 
tary sequence pairs over other forms of noise modulation, Figure 1, 5-g 


shows the same amount of delay as did all the other forms of matched 
filter detectors, and shows that the sum of all the energy sub-pulses is 
incorporated into the main pulse at exact match, 

Each of the systems described in this introduction has had the char- 
acteristic that the transmitted signal is spread over a wide frequency 
spectrum, This has in itself two advantages: One is that enemy detection 
of the signal becomes more difficult, and the other is that even after 
detection the jamming of a wide spectrum noise-like signal becomes 
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extremely difficult. Since the characteristics just described are de- 
sirable for a secure communication system as well as fora radar 
system, all of the advantages that accrue for the complementary se- 
quence modulation in radar would apply in a communications scheme, 

The initial major objectives of this investigation were two in 
number, 

Akg The formulation of a set of operations on complementary 
sequences which form a group, along with a theoretical 
study of the invariant properties of sucha group, 

2, The search for new complementary codes, An exhaustive 
computer search for a new code of length 26 was made 
as wellasa partial search for codes of length 34. 

During the investigation two more major objectives were added 
to the work, 
These were: 
5 The application of Hamming distance to both the comple- 


mentary sequences and their group formulation, 


4. The application of Hamming vectors to complementary 
sequence pairs and also to their group formulation, 
The original objectives were oriented toward the generation of 
new codes, while the last two objectives were oriented toward the 
decomposition of codes. All of these objectives were attained and 
consitit ute a large portion of the rest of this work, The next chapter 
is a study of complementary sequences designed to bring the reader 
to the state of the art, and contains some proofs which have not ap- 


peared elsewhere, P2595 Ae 


Chapter 3 is the formulation of the 
Operations group on complementary sequences and, along with several 


theorems and proofs it also contains the identification of this operations 
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group with a known group, Chapter 3 as written requires only a min- 
imum of group theory knowledge. Chapters 4 and 5 have to do with 
the Hamming distance and Hamming vectors of complementary sequen- 
ces both in standard form and in group form. Chapter 6 isa short 
chapter and is concerned with the proof of two theorems used in the 
searching for new codes, Chapters 7 and 8 describe the computer 
search for new kernels of length 26 and 34 respectively, the actual 
computer programs used being listed in the appendicies, The last 
chapter, Number 9, contains the summary, conclusions, and suggestions 


for future research in this interesting field. 
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CHAP bh Roi 
COMPLEMENTARY SEQUENCES 


This chapter is presented to provide a background in complemen- 
tary sequences to enable the reader to better understand the material 
in later chapters, Most of the materia] in this chapter is an adaptation 
of Golay's ''Complementary Series'' with expansions and deletions to 
fit the needs of this paper. ? The formalization of some proofs and the 
addition of new proofs given here have not appeared in print elsewhere; 
however, it is obvious from the tone of Golay's work that he was aware 
of these proofs, This modification of Golay's work is offered for com- 
pleteness of the dissertation and not as new work, 

In the introduction a pair of codes of length n were discussed which 
had the characteristic that when detected by matched filters, the sum of 
the two filter outputs was everywhere zero, except where the signal 
patterns had zero time delay with reference to the filter patterns, At 
the zero pattern delay time the output was 2n, leading to an infinite 
coding signal-to-noise ratio, <A study of the necessary and sufficient 
conditions for this to be true will now be considered, 

Assume a complementary code pair (A,B) and suppose that one of 
the pair, i.e. A, is longer than the other. If this be true then there will 
be an output from the A filter with no cancelling output from the B 
filter, This is not allowed; therefore the A and B codes must be of 
equal length. Therefore assume a code A, n bits in length, with 
elements a, through an where each a. is either +l or -1l, Also, 
assume a second code B, n bits in length, with elements b. similarly 
defined, Detect each of these codes with a matched filter (the time 
inverse of the code) and then sum the outputs of the A filter and the 
B filter. Attime 1 (see Figure 2,1) the signals a; and b, have 


entered the filter segments a and De respectively, Since at this 


)3 








time 
be care b3b,b, 
The a. and b. of the filter are 
either + 1 indicating straientetoes 
amplifiers or inverters, 
fare») A 423458) signal 
a Ao] Ani2eeees 2 Fy filter 
b3b,b, signal 
by bao Da ee es -b) filter 
iaraxve 2 a 4a,a5a) 
aa ane | ens ake neler 
b,b2b5)) 
De ee l ae Z 
Mire: I A. sé ce os a 
r 1 
a 


iui cela Ant]-(r41) Ft 


Matched filter detection of a complementary 
pair of length n, 


PIGURE@Z- 


14 





time the delay between signal pattern and filter pattern is not zero, the 


summed output voltage must be zero, This gives the condition aja t+ 


b)b. = 0. Attime 2 (Figure 2.1) the signals are one segment farther 
into the filter and the output must again be zero, giving aja,., + 
a,a_ +t b Dre ti bob. = 0, 
Ewteetime © 3: a) ay gtasa,_)ta3a,+b)b _5+b5b__ ;tb3b_ = 0 
At time r (r # n): ate 2 eee fete Tad toe tees 
ron les Peo 
+ bb 20 
eat ot 
mat time n: a,ay ta,a,+......+ ann +b) b, eee + bbs = 2n 


for at time n the pattern of the signal must exactly match the pattern 


of the filter. The explicit statement of the conditions given above for all 


” 
= 


times is 
i=n-j i=n-j 
fe a aiainy F am ree ae Lie OMT ps. 3) 
— 7 =a j -= 0 


An examination of the equation at time 1 (j= n- 1) shows the following 


eight possible solutions: 


ay ay by a. 
l ] l - 1 
] ] =o 
- 1 -1 ] -1 
- 1 -] -] J 
- 1 iE ] ] 
] -! i | 
- 1 l =" = 1 
] -] = | - 1 


These possibilities show that if the A pair are alike, the B pair must 
be unlike, andif the A pair are unlike, the B pair must be alike, The 
equation at time 2(j=n-2) indicates that if the number of likes (like pairs, 
i.e, a, and ay?) of the A code and the number of unlikes (unlike 


pairs) of the B code are equal, andif the number of unlikes of A 


equals the number of likes of B, the equation is satisfied, Extending 
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this reasoning to time r shows again that if the likes of A equal the 
uhlikes of B and if the unlikes of A equal the likes of B the equation 
is again satisfied, 

It will now be proved that the last half of the previous statement is 
not necessary, since the first half of the statement contains the necessary 
and sufficient requirement. 


Let UL = number of unlike pairs in A for the spacing specified. 
(j in equation 2, 1) 


Let L, = number of like pairsin A for the spacing specified, 
(j in equation 2, 1) 


Let UL = number of unlike pairs in B for the spacing specified. 
(j in equation 2, 1) 


Let L, = number of like pairs in B for the spacing specified. 
(j} in equation 2, 1) 


Assume that L, = UL: we want to prove that oa = ine 


Since the sum of the number of like pairs at each spacing with the num- 
ber of unlike pairs at the same Spacing must equal the total number of 
possibilities, 
U_ +L, =n-(n-r) for r #n where n is the length of the code 
a a . 
and (n-r) is the number of units of time delay 
from exact match, 
For the same reason 
UL + L. =n-(n-r), 
But since UL =L~, 
a 
C.. + UL = n-(n-r). 
Therefore U_=L. 
a b 
The definition of a complementary pair of sequences will now be 
given, A pair of binary sequences of equal length with the number of 


like pairs of one sequence equal to the number of unlike pairs of the 


other sequence for each possible spacing is said to be a complementary 


code, 
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Two schemes of representation will be usedin this paper, These 
schemes will be used interchangeably at the convenience of the author, 
The first representation is the one already presented where each a. 
and b, is either +! or -1, The operation used in this case is 
ordinary multiplication, The second representation will use the opera- 
tion of modulo two addition denoted by 8, where the elements a. and 
b. are either O or 1, That these two groups (1,-1;¢ ) and (0,1; ® ) 


are insomorphic is shown in the following tables. 





g 0 l 

0 0 ] 

l | 0 
lL = 0 
= = } 


The first form will normally be used when correlating of the codes 
with their matched filter is under discussion while the second form will 
normally be used when discussing the intrinsic properties of the indi- 


vidual codes. 
The necessary and sufficient condition for a pair of sequences to be 


complementary, is interms of the modulo two representation 
i=n-j i=n-j 
F.c- alwasw..) = 2 (b.9b., .91) for all j, ]<j<n-1 Ze 
p22, iMag) = Z| (oieb yy OY =) (2, 2) 


or an equivalent formulation by change of summation index 
1=j 1=j 


y 
os (a.® a = =] (b.@b 


j i=l *°i n-jti 8 1) (2, 3) 


n- jti 


The above % is standard summation, If modulo two summation is 
desired the symbol ¥ will be utilized. (The equivalent conditions are 
given since they are both used in formal proofs ina later portion of 


the paper, Due to the possibility of confusion it seems more appropriate 
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to introduce the second form now rather than later when needed, ) 


oe = L, has been used to represent these same equations, and will be 
used elsewhere in the paper when it 1s more convenient to use it, 

An interesting property of any system of codes is the possible in- 
variances it may have under various types of transformations, It 
therefore seems appropriate to consider the possible transformations 
ona complementary sequence pair which leave the pair complementary 


and of the same length after the transformation. A necessary tool for 


this study is the property that a ai, = ake a. 


: where the symbol 


+3 
a. means the complement of a.» O=1, l=0, 
Theorem 2, 1 
The modulo two sum of a pair of binary numbers is equal to 
the modulo two sum of their complements. 


The truth of the theorem is obvious from the following exhaustive 


table of four possibilities: 


a. § a a. & aa 
1 ® 1 = 0O ae 1 es 
1 ® O = 1 1 6 0 S=550 SO 1 = 
O we YT = 1 0® 1 = 1 @¢ 0 «¢= 1 
0 @ 0 = O 0@ 0 = 14 1 = O 
Therefore a. Gay; = a. v ij , 
2K aie 2 


In each of the following proofs it will be considered that A and B 
are acomplementary sequence pair of length n. Complementing each 
element in a code is called complementing the code, 

Theorem 2, 2 
Complementing the A code, or the B code, or both the 
A and B codes, results ina pair of complementary codes, (The 


proof is given just for the A code but is identical for the B code and 
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from both of these for the A and B codes.) 


By hypothesis 


1=n-J i=n-) 
F.= 5 (a@a,) = DZ (b @bi,, 01), 
i= | c—" 


Complementing the A code transforms 


aeaatt a. , 
i i 


but by Theorem 2, l 


So Rare = oe : 


Upon substitution 


i=n -] 1=n-] 
=e Dy 5 5 om Z : 
e - (a. Y ae - (b. 8 Bont ® 1), 


Therefore complementing the A code has no effect on complementarity, 
and similarly for the B code ~ for both codes. 
3 SS x 
Interchanging the first and last bits of a code, the second and next 

to last bits, the third and third from last bits, and so on is called time 
inversing or time reversing a code, 

iieoreimc. > 

Time inversing the A code, or the B code, or boththe A 
and B codes, results ina complementary pair. (This proof is again 
just for the A code, but also applies to the B code and therefore to 
both the A and B codes, ) 


By hypothesis 


= a 2 | 
1 (a. a, a (b @b.,. @ 1). 


.) 
J h—oE ¥] T= 1 J 


The time inverse of A causes each 
a.r7ra i. 
iL n+l -i 


Applying this time reversal just to the A code shows that 
i 





l=n-]j 1=n-]j 


a 7 (2, TA - ae (aaa *nt1-(itj) ) 


jlo 


Expanding both expressions for a few terms and remembering that 
modulo 2 addition is commutative as is regular addition we see that 
(on the left) (on the right) 


for j=n-l, F174) 84, Fo} = 3, % 4) 


ror J=n-<; Pa-2 7 219 ey + a4® a, F277 4,8az,ta Ga 


n-l | 


These expressions indicate that by a change in summation index the two 
forms of es are equal, as is the case. Therefore time reversing the 
A code has no effect on the complementary property, and similarly for 


the B code and for boththe A and B codes, 

The altering of a code is a transformation where every other bit of 
both codes is complemented. Ina different portion of the paper a dis- 
tinction is made between altering odd bits and altering even bits. In the 
following proof, however, just the even bits are altered, although a 
slightly modified proof would hold for altering the odd bits. 

Theorem 2, 4 
The result of altering a complementary pair of sequences is 


again a complementary sequence pair, 


Since the codes are a complementary pair, 
i=n-]j i=n-]j 
a = 2 2 ta 
F > (a. Ba. (b. oo 1. 


P) 
el T) al 


Divide the ig into two groups, FL for j= an even number and F 


for j= anodd number, The portion of the above equations for j= an 


even number is 
icone fan= Kk 
=e: 
a a a! | 
i=] t=) 


® 1). 


0) 





For the theorem to be true, when j is an even number the above set 

of equations FY must be satisfied for all k when the pair is altered, 
(The proof is given only in terms of A, but duplicate steps must also 

be applied to B.) 

When i=an even number, then itk=an even number, In the alter- 
ing process since all even bits are complemented a. ® aia, 7 a. Oa.) ; 
but by Theorem 2,1 these are equal, When i=an odd number, then 
itk = an odd number and these are not affected by the altering process, 
since only even bits are complemented, Therefore F,. is invariant 


k 


under altering, or symbolically F. = Fait’ 


Examination of F where r is an odd value of j, shows that of 


the pair as a. one of the bits is complemented and the other is un- 


ra 
changed. 
7 : Gite (a,® Site OF She a ity! a : Sa aan 
Therefore ae a ee (n-1)-F for both the A code andthe B code. 


This signifies the likes of the A code for odd spacings are changed into 
unlikes, and the unlikes of the A code for odd spacings are changed into 
likes; however, for each change in A there is an opposite change in 
B so the total likes in A _ still equals the total unlikes in B for each 
spacing, It is to be noted that altering both the even and the odd bits is 
the same as complementing both codes, 
oe oe XK 

The operation of interchanging the A and B codes, although 
trivial, is also a transformation, There are a total of sixty four pos- 
sible transformations, These can all be generated by combining the 
above listed transformations, They will be discussed much more thor- 


oughly in Chapter 3, 


il 





Whenever a new facet of a pair of sequences such as the comple- 
mentary property is discovered, it is natural to wonder just how general 
this property might be. One question of interest would be the possible 
length limitations of the code pairs that might be forced by the comple- 
mentary property. The following proof due to Golay shows that 
complementary sequences must have a length n of the form eater 
where n is an eveninteger and x and y areintegers, The first 
theorem to be proved demonstrates that n must be an even integer 
and the second theorem shows the elas Low aae 

Theorem 2,5 
A necessary condition for a sequence pair to be complementary 
is that their length be an even number, 


By hypothesis 


i=n-] i=n-]j 
p= oy aQa... = od b.@b.,. @ |} ’ 
iz] 7? 7*J i=l) ae 
leet. = F. modulo 2Z, so that 
J i=n-} 1=n-]j 
C. = ) a.@ ee = E b. eT ee lo 


i=l 


An expansion of G fora few different values of } will suggest the 


proof, Start with j=n-l. 


Gi] = a 9a, = b) 6b 8 Le 
G> = a Ga _ 042%, = b eb. | 91 8b,%b el, 
G3 =a Ga __ 78a, 1 92, - 1 93-394 =P) Ob, 28 Ob. | 218b2 9b, 61, 
Change the order of the sums in G2? and G3 to achieve the following: 
2 
G24] ®a,ba_ 82,2) 9b, eb, 9b, io) 5 ae 
3 
G 374 | 947 9a,8a_ 78a | 8a ,=b Ob, 9b20b 4b | 2b, z [te 


In general for j=n-r 


G7 98 29. = Oa Oa (rt) Oa (r42)" > Ga. | Ga ,=b) Ob, 9. gate 


(aya 





c 
62 OOD, (4419, (742) + OP 19,9 EL, 


ancuon aaj—m=(r+)) 


=a,9a,9...8a 8a Ga 6...@%a 
1 Zz 1 IS 


ne er ores Ga =b eb, 6. ie 


rt+l 
ae Ob Ob. 0b, _ 9d, _ (eee omen »o @ ¥ 


ert) n-] 


l © 


Forming the sum GG. ~? G will leave only five terms in the 


n-(rt+1) 


equation since G+ and iret) differ only in Aner? jel]? DEL? 


and a, bcause 


lu ct | 
o 10 y l=l, 
] ] 
n-r easier ta) ~ eeleaner: ~ ree: Onn L. : 
Adding B41 85L_+ modulo 2 to both sides of the equation gives 
Bee etl nes ‘I. (Gs 


Let n=2s-1 and r=s-1; thento satisfy the equation just derived 


a ®a @b . Ob 


s 2s-l-stl = a 8a 0b Ob .= 1, which 


2s-l-stl 
is obviously false, 
Therefore n cannot have the form 2s-1 and must be an even number. 
x 1. a5 
Equation 2,4 which was a step in the previous proof, is a very 

important necessary condition for a pair of codes to be complementary 
and will be referred to in the future as the parity check. The next 
theorem which shows the nee form, starts with the basic assump- 
tion that A and B are a complementary pair of length n. 

Theorem 2,6 

A necessary condition for a pair of codes to be complementary 


is that their length be the sum of the squares of two integers, 


By hypothesis 


i=n-j l=n-j 
a = fa (a. are = a 122 ng! Z 


Assume that the A code satisfies a and that the A sequence contains 


(a) 





p ones and n-p zeros, The A code has aweight of p. Similarly 
assume that the B code satisfies ee and that its sequence contains 

q ones and n-q zeros, The B code has a weight of q. For the total 
likes inthe A code, each bit is matched against each other bit exactly 
once, For the total unlikes in the B code each bit is again matched 
against each other bit once, There are two possibilities of likes in the 
A code, 1 matched against 1, and 0 matched against 0. The total 
number of like pairs for the A code is therefore the number of com- 
binations of p ones taken 2 ata time plus the number of combinations 
of n-p zeros taken 2 ata time, for the unlike pairs inthe B code, 

q ones are matched against n-q zeros. 


2 


This gives the equation 


pip-1) , (n-p)(n-p-1) — 
oe 2 


expanding and simplifying 


q(n-q) 


penceentea- + 2p--2np, 
adding and subtracting 2pq gives 

er eames ten 2 enpt2qp-2qp, 
collecting terms gives 

m= A areae +2qp-¢enp-2nq + pe -2pq hee 
combining terms gives 

n= (n-p-q) + (p-a)° (2.5) 
Therefore complementary sequence lengths are permissible only in 
lengths which are even numbers and formed by the sum of two squared 
integers, 

ok pS ok 


A list of all possible code lengths upto 200 with the number of 


ones allowed in the A and B codes is givenin Table 2,1, 
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*Indicates possible kernel 


Code 
Lengths 


«2 
4 
8 

*10 
16 

*18 
ZU) 

*26 


Ore. 
*34 


36 
40 


a) 


Be 


#5 8 


64 
68 


‘2 
*74 


80 


oc 


*Q 0 


*Q 8 
100 


Unordered 
Weights of (A, B) 
(2,1)(1, 0) 


(3, BONS Ue 
(6, 4)(4, 2) 

(7,6)(4, 3)(7,4)(6, 3) 
(10, 10)(6, 6)(10, 6) 
(12, 9){9, 6) 

mee 11) (9; 7)(141,,.7) 


(16, 15)(11, 10) 
(15, 10)(16, 11) 


(20, 16)(16, 12) 


(21, 18)(16, 13) 
feito (is. ro) 


(Zr 15)(21, 15) 


(24, 22)(18, 16) 
(24, 18)(22, 16) 


(30, 25)(25, 20)(29, 28) 
(29, 22)(28, 21)(22, 21) 


(31, 27)(31, 25) 
(eee (27, 21) 


(ea e)(27,24) 
(ae 72. 24) 


(36, 28)(36, 36)(28, 28) 


(39, 37)(31, 29) 
(Some) (67, 29) 


(42, 36)(36, 30) 


(43, 38)(36, 31) 
(43, 36)(38, 31) 


(46, 42)(38, 34) 
(46, 38)(42, 34) 


(46, 45)(37, 36) 
(46, 37)(45, 36) 


(51, 48)((51, 42) 
(42, 39)(48, 39) 


(56, 49)(49, 42) 


(55, 55)(55, 45)(45, 45) 
(57, 51)(57,49)(51, 43) 


(49, 43) 


Code Unordered 
Lengths Weights of (A,B) 
104 (58, 56)(58, 48)(48, 46)(56, 46) 

*106 (60, 55)(51, 46)(60, 51)(55, 46) 
116 (65, 61)(55,51)(65, 55)(61, 51) 

*122 (67, 66)(56, 55)(67, 56)(66, 55) 
128 (72, 64)(64, 56) 

*130 (72, 69)(72, 61)(61, 58)(64, 57) 

(73, 66)(73, 64)(64, 57)(66, 57) 
1 36 (76, 70)(66, 60)(76, 66)( 70, 60) 
144 (78, 78)(66, 66)(78, 66) _ 

*146 (81, 76)(70, 65)(81, 70)( 76, 65) 
148 (81, 79)(69, 67)(81, 69)( 79, 67) 
160 (88, 84)(76, 72)(88, 76)(84, 72) 

* 162 (90, 81)(81, 72) 

1 64 (91, 83)(81, 73)(91,81)(83, 73) 

*1 70 (94, 87)(83, 76)(94, 83)(87, 76) 

(92, 91)(79, 78)(92, 79)( 9 er) 

* 178 (97, 94)(97, 84)(84, 81)(94, 81) 
180 (99, 93)(87, 81)(99, 87)(93, 81) 

*194 (106, 101)(93, 88) 

(106, 93)(101, 88) 
196 (105,105)(91,91)(105,91) 
200 (110,100)(108, 106)(100, 90) 


(94, 92)(94, 108)(106, 92) 


All possible complementary 
code lengths up to 200 with 
unordered possible weights, 


TABLE 2,1 


Vad, 





The shortest possible complementary pairis A=l1, B=!0, This pair or 
any of its transformations is called a kernel of length two, or the quad, 
A kernel is a basic length code which cannot be decomposed into shorter 
length codes by an inversion of the standard generating methods to be 
explained later. Some possible kernel lengths are 2, 10, 18, 26, 34, 
50, etc. although codes for all of these do not exist. Complementary 
sequences which are not kernels are called composite complementary 
sequences, 

The above list of kernel lengths did not include n=4 or n= 8 
which are possible complementary sequence lengths, These if they 
exist must therefore be composite. All possible codes of length four 
can be generated from the exhaustive list of possibilities for four binary 


digits given in Table 2, 2, 


0000 1000 
0001 1001 
0010 1010 
OO VST Th 
0100 1100 
0101 1101 
0110 1110 
O1ll ele) 


Table 2.2, All Possible Binary Numbers of Length Four. 

The unordered possible pairs of ones in the codes are (1,1), (3,3) 
and (3,1). Thus all possible codes of length four have either one 1! or 
three 1's because of this limitation, This reduces the table to the 


eight following numbers of length four: 


0001 Olli 
0010 1101 
0100 1011 
1000 1110 


The possible unordered (1,1) code pairs are the following four: 
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The possible unordered (3,1) pairs of codes are the following eight: 


O11] Oa Tbe 1110 1011 1011 1101 1101 
0010 0100 0100 0010 0001 1000 0001 1000 


The possible unordered (3,3) code pairs are 

0111 0111 1110 LO 

1011 oT POLL Prod 
Comparing the last code to the kernel of length two, A=11, B=10, 
shows that it might have been constructed by writing in time sequence 


S; = AB=a,a,b)b, 2 ILO 


= AB=a.a.b.b 


So i222 = 


EO 
The next to last sequence pair might have been formed by inter- 
lacing A and B inthe following manner: 


= a,b a,b, =F ele O 


T3 = a)b,a,b, = 


All the rest of the 16 pairs could be considered either as transforma- 


Or 


tions of the last pair, or as being formed by the same two composite 
generating operations used above on the 8 transformations of the ker- 
nel of length two. A general proof wili now be given to show that both 
the time sequence scheme (S,, Se) and the interlace scheme (Tj, Ts) 
will always form complementary sequence pairs providing that A and 
B are acomplementary sequence pair, 

Mineorerm 2. 1 


If A=a)a542. oe ee 


B=b)b,b,. ecdenns slo, b 


are a complementary sequence pair, then 


On Sy a Le b 


n 


v= aaa. ; ab) b,. en be 


are a complementary sequence pair, 


at 





Theorem 2,8 
if (A, B) are a complementary sequence pair, then 


@e= abn ae bays ss ie a_b 


imglin 2: 2 nn 


D=a 
form a complementary sequence pair, 


The following notation will be used in the proof: 


1=n-j 
ie => a.Oa..,. for-alle;) < | < 1am 
a izl 1 it] 
_ ifn-j : 
L, = 5 Sa aie for all ja ere 
i) 


It is understood that the two CD sequences in the two theorems 
above are not the same but because the proofs are identical two ete of 
symbols will not be used. For.each spacing j there are three possible 
ways for bits to match, within A, within B, or froma bit of A toa 
bit in B. These possibilities are denoted by L. Ly and Lopor U_> 
UL, ee respectively. 

The necessary and sufficient condition for C and D tobea 
complementary sequence pair is given by L. = L,+U +L, and 


U.=U +U- +t U - where L =U 
a a Cc 


d b b d° 
a aa ~ NE gnc Li, ~ ae by definition of complementary. 
Adding 
L, + Ly = U,tU,» 
cae i=n-l i=n-1 Z 
UL = UF since ee eee = Ss b 9 bit; 


so substituting 


b 
All bit matches from A to B inthe C sequence are paired inthe D 


ee teem Oa t UE 


sequence with bit matches from A to B therefore, 
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L., =U_>- 


ab ab ’ 
and adding 
Pee, Fe UL eet UL or L. =U, 
at xe ss 


Golay describes two other methods of generating composite sequen- 
ces from shorter complementary sequences. These methods will not be 
used too extensively in the rest of the paper and will therefore just be 
mentioned rather than proved, Given two complementary sequences 


pairs (A,B) (C,D), A of length m and C of length n, 


Let 
U, = ASIAS? .., ...Acn Bt Be Ban 
Us = A™,....00.....A7 “BN ,,..........B%, “ 
and 
eo dD vec ascastas Shem Bon 
PA Bo? |. ook... cease es ee eee es ah 


It can be shown that the pairs of U are complementary as are the pairs 
of V. Where, if an exponent is aonethe A or B code is left un- 
changed, and if the exponent is a zero the A or B code is comple- 
mented, The lengths of codes which can be generated from these two 
methods are 2mn where m is the length of the (A,B) sequence pair 
and n is the length of the (C,D) sequence pair. There is also another 
special method of generation which applies only to codes of length lan 
This method is explained in Chapter 5 where it is necessary for the 
completeness of a proof, 

The kernel of length 2, or the quad, has been discussed in the 
previous pages. There are two kernels of length 10, These are 
1001010001 and 1000000110, anda second Kernel 0101000011 and 
0000100110, The next possible kernel size is 18 and Golay has proved 
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by exhaustive search that no kernel of this length exists, Kruskal has 
since then completed the proof analytically. a One of the goals of the 
writer in his study was to make an exhaustive search for codes of kernel 
length 26, This search disclosed there was only one kernel of length 
26, Chapter 7 of this paper describes this search in detail and 


Chapter 8 describes attempts to find kernels of length 34, 
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CEAPIT ER il 
Poreorwour OF OPERA Pion. 


A large portion of this investigation is devoted to the set of allow- 
able transformations on complementary sequences, These allowable 
transformations possess the characteristic that the length of the sequen- 
ces and the complementary property are invariant under transformation. 
To be useful in searching for new kernels it is also desirable that these 
transformations have the following features: 

1, That they forma group for all general transformations, 


Where a general transformation is one which can be 
applied to any complementary code pair. 


2. The transformations as represented be simple manipula- 
tions on the elements of the group. 


3. The supplementary characteristic (to be explained in 
Chapter 6) be a property of the representation. 


4. The parity check be easily made. 

Some of these transformation operations were discussed, and the 

proof of invariance given for all group generators, in the background 
chapter, Chapter 2, on complementary sequences, As a continuation, 
a brief review of the transformation operations anda tagging of symbols 
to these operations will prove useful. Each basic operation symbol and 
its definition is given in Table 3,1. More definitions are given than are 
necessary to generate the group; however, the redundance seems to the 
writer to be an aid to understanding. 

A check through Table 3.1 shows that the transformations affect 
the complementary sequence pair in four different ways, C,;, Ce, 

Ti, Ts affect all the bits of one of the code pair, E,T,C affect all bits 
of both codes, while A; and Ag affect half of the bits of both codes, The 
identity, I, of course has no affect on either code, A reasonable way 

to form the elements of the group to account for the above affects would 
be to divide each code of the pair into two units, These would be A, 


(even bits) and A, (odd bits), also B, (even bits) and B, (odd bits). 
Sil 





Th 


Te 


Az 


Ci 


C2 


The Identity 


A code time reversal 


0 ee : 
1 nt+l-i 


|Saeiaad Sy 
1 


1 


B code time reversal 


Deas as 


1 n+1-1 


Time reversal of both codes 


ana : 
il nt+1-i 


Di ee 


Alter odd bits 


See el 2S SD. 


ee a 
Alter even bits 


tOi-1 “eae Zien 


bee bea 


b Daan ey nes 


Z2i- 1 
Complement the A code 


Complement the B code 


Complement both codes 


oe 
Some operations, 


as 
Exchange the A code with the B code 


a. 7 b, TABLE 3,1 
b. yy S- a2 





il 
© 
ey 
© 


OD WB YP Pp 
1 iT iT 
oT Pp 
tN 
oT 
ds 
oT 
= 
iy) 


= bobyby. wee ce ee bd 
However, this particular division for representing the complemen- 
tary sequences proves to be very awkward for the time reversal opera- 
tions (T,,T2, T) and for the ease of the parity check condition, 
Fortunately, one simple change, that of time reversing the second and 
fourth row (even bits of A and B), clears up these difficulties, a 


The method for converting from a complementary sequence pair 


into the standard group form is therefore 


ha a ee OO el 
= eae ee rr a 
hoe 3 n-l n> iz 
= ha enone (3.15) 


Li bee 


B=b, bee... ...b bX 
Laz. 3 n-lon v=o ae 
n n=2 Zc 
If the matrix (I II III IV) is taken as the symbol for the above grouping 
the Identity operation, I, should yield (I II III IV). Post multiplying 


meet if] IV) 


by gives (I II III IV), therefore I can be taken 


Qr Ga 
oO Oo =o 
oa GS 
=o Ceo 


as the identity matrix of rank 4, using the operation matrix multipli- 
cation on the right. Similarly T; has the effect of exchanging the first 
and second columns yielding (III III IV). Now post multiplying 


(I II Ili IV) by 


Cale O° 60. 
; ; : yields (II I II IV) the desired result so the 
oO 0 | 


matrix just used could represent T,. Similarly, 


a2 





1 0 “Oe 0 1 0 0 iid 0 1 0 0 0 
eo ter 0 «+. 1 0 0 0 ~-2. Oto 0-20 1 0 0 
aCe 1 CC ee 3 oo 8 20 8 Co 

0 0 1 0 0 0 1 0 0 0 0 1 0 0 0 1 

1° @ “te 0 oO 0 Oe 

_ (io 0 0s ecm) ee py oe 
caeeeceo 1 0 “= 9 60 10 @* 4 6 Deo 

00 0 1 0 0 0 1 Oy sles 2A 


Matricies are often a very convenient method of expressing trans- 
formations. However, when the elements are few and the numbers in 
them simple it is often more convenient to express the operation dif- 
ferently. For instance (I II III IV) operated upon by Ty yields 


(II I Ill IV). This could be just as well expressed as 


Sele it 1V 
ele tlh 1V 


at 

This indicates that the top row is transformed into the bottom row by the 
Operation, or in even more shortened form, Ti = (II I III IV) with the 
top row understood which means, as is common practice, 1T, (I II III IV) 

= (II I Ill IV), Another example is A, = (I II Ill IV), [which means 
A, (IIL III Il) = (1 II ll IV)] where the bar across I and II indicates 
the negative element in the matrix or in actuallity the complement of the 
portion of the code contained in I and II, Table 3,2 is a list of all 
32 transforms generated from Table 3,1 both by symbol and group 
element representation, There are actually 64 transform operations 
in the entire group. However, the exchange operation (III IV I II) is 
not considered along with all the 32 elements it would generate, In 
this paper the pair (I II) and (III IV) are generally considered as 
being an unordered pair, which eliminates the exchange operation, 


The second feature listed under desirable characteristics, that 


of simple manipulation of the elements of the group is satisfied. 


The supplementary property which is discussed in Chapter 6 is 


also satisfied by the formulation of groups from complementary 
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sequences by this method, The parity check, a %b. a ®b 


ales nels ee 


is also satisfied under all transforms and is easily checked since each 


column of equation 3.1 is one solution of the parity check. 


1 n- | 

I= aan .2 o+++ 8nd] 22 (3.1) 

I= BBD. os +>. b.....b | 

tate - 2° nti 12 
ieaee{1 1 Il Iv) R= 4! Il tv 
Ai= (I Il Ti Iv) @r= 1 ray 2 a 
fe= 1 Il Ill IV) P= (De 1 ei 
ie= (il 1f Ill IV) O= (ies lie 
To= (1 II IV III) N = (i pes) 
toot. 1 IV III} M =e Gy enn) 
Ga iil iil Iv) ; Le = (10, Tee) 
ee-nl Il Ii IV) K= (1 ipay 
ero) DL Ill Iv) J= (1 I IV Wy 
Z= (II I Ill Iv) H = (I II IV Ty 

Sit Tit IV) G= (1 I Iv Dy 
= (II I Ill Tv) F = (I Ir Ii Iv) 

W = (II I III IV) D= (I II Il IV) 
Ve= (1 Il Iv Ip G@= (Il 1 ty 
Weil tl iV 11) B= (1 I IV Ti) 
S= (I II [V IID mn = (IIT IV IQ) 


E =(III IV I I) 


Elements of the unordered operations group. 


TABLE 3.2 


35 





PGE 


62) 
iv ay] a et — 
OCH HOMEDZORSEDNOHKENHAMHYIHOOURUaLZaH 


CQ aoiec -t 
aero Sse OS SSO Nes HOON eee Hee eo 
MOM PRORBHEMMHEAZHONOOPNND SAM aH HONEA 

oa YN 4 qo \ 
AxttHOUFD FL4ZAOCOHDMHAMOAUOUNMN NH EMF AOHH HO 
ay) os = 
ft COQOOHMMMOBBZIS MMH OOANKME ED ENQHOHKASA 
eo =, 7 
OHHRAENZEOODNES EN HMOAMAMEm Be OGHUMA aD 
ay) iq] ot —t 
TM ASS OAOBZOMA ETA AEN KH HOOOOH QB DME 
ay) ca ot 
DOOHAUMMAOPRMNDNS <Ha eee eH DR OUMERAD 
a; ay) o 4 
MOS HRD HOAMOQEMH tO tH Ss OAOOHHEOBZENANKO 
rc ay] ra YQ ea ( 
ABZORFHOHMAKE aS EOMADRPOORHHHUEKNOES 
ay] rl ay) — 
SeMON>bOOOMEOAMAMSARNHHHOeUZONRDREHA 
al a) M ral a ray) 
ZAAOMNOOAMAAAHRAMRENEOHHOOS HM EHWDeE 
4 qd a © at] O 
OQOHRPWMAN MOOS ART OOFA MEP NHOUBROUOQMHA SH =ebpk 
~O mM ~ ay} 
MEMOHHNK AHS MME SIOHOOD Re SENaHOoOROOE 
ay] ; = 
eee ER HEN AMES MOND SZOH RD HNM ae HaAuoon 
ra] — +t 
Bre NOME OME HHAHOODHO IRON EAE ANDES Zore 
rat] ay} a 4 
NMDoDSe> StU KHORAHOOOHOZEMAQKNE SE aME YH 
N a) aoe a) 4 
MMW>meaeM FO HOM OAHOUDUUEBWAAOt Mea at eM STrOseN 
a a) ct gy) ay) co 
Pa DAMS KtnNnNanroktOnOUBZOMERaA aA KOHN AaVWes 
ay) Q ea 
eS MN HOBO EDOOHOHOHOEMHASARMBRANnZOmsS 
a) 4] = oA ay] 
MEHMNZOMNNOOOHOAHOH MS OHREDOPMODYQE Mee 
= ay) ~ a MN 
MN M EO eM RPHOOOHAH OHM AAOQEDAS OSSIAN 
ra) ay) ot ay) 
NH EHMOBZ ZOE OHOOOHOEHSDHOAMAANEREAHMD 
MQ et ca CQ 
HOOK MOHEHNDE Re HMHNSEe ee Se TOMYQOBZOEMAMD 
M oi - a) rea] a) 
ae ase meee Ea moma eet Se HOMO O Me Ce eee] 
ay] a) + | 
OBS ON EHH em eOAMAEORDOOUKUOH MOESNME NE 
iq) fay) OQ a 
AOmMAOHMNEARZEOMS HE MDH HE SHDONE Or EOORUm 
fay} a | 
AmODHONEQHOUZHEM SHH HEE KNADMOOUURE 
rc — 
OG GHA EOMHMENHNEEPDMUOZOMAHH SOAR ENMY 
ay) NQ eae 
SOMO Ome MOE KHN DEON EMQNEMOB OHH ae oe 


o ay) rq) 
SOHOOMAOSOHNEMBENDE MERRY Z MEME aaah 
ray; 


oA NU A 


TABLE 3.3 
Operations Multiplication Table 


fay) 
FOOOdt eR HHN HM eb DpDnmomoFr a oIM> Hoo eam eae 
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Operation 


— 
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(Side) (Top) 





Another feature inherent in this formulation of the group, which is 
very useful in searching for new kernels, is the interlace breakdown 
into the group, This is useful because a kernel may be considered as the 
interlace of two half length codes where, as will be shown in Chapter 7, 
the number of zeros and ones in each half length code are easily pre- 
determined, 

Although this group of operations is indeed a group, it has never 
been proved formally. In order for a collection of elements to be a 
group, they must satisfy the following four S oaclineines- 


. <A group has closure, 


i 
2. A group is associative under the group operation, 
3. A group has an identity element, 
4. Each element in a, group has an inverse, 
That this operation group has closure is demonstrated by Table 3.3 
which is the full multiplication table of this group of order 32, The 
elements under the operation multiplication are associative since they 
are expressible as matricies and matrix multiplication is associative, 
The element I is the identity, and as seenfrom Table 3.3 all elements 
have an inverse, Therefore these 32 elements have satisfied all the 
requirements for a group, 

There are three special subgroups of order 4 which within the 
subgroup deal with only one type of operation, T, T2 TI, A; Aza Cl, 
Ci C2CI. All of these subgroups have the unusual property that each 
element in the subgroup is its own inverse, This group of order 32 
has been identified as isomorphic with Senior's group number 44, te 
This will be discussed in some detail at the end of this chapter. 


A study of complementary codes in (I II III IV) form, or as we 


will call it sequence quadruple form, reveals many interesting proper- 


ties, Some of these properties and their formal proofs will follow, 


a 





In these proofs the complementary pair (A,B) will be assumed to be of 


1=j 
length n, The symbol L, will mean a a taney i @1) for all 
‘ : : : _ 1s 
j, 1S js n-l and similarly Un > hoe 
a 1 nti-j 


at 
Theorem 3,1 


If (A,B) are a complementary sequence pair and if either (I,II) 
or (IIJ,IV) in sequence quadruple form are complementary, then the 
other pair must also be complementary. 


1. Assume (J,IIJ) are comlementary. 

Gis L. = Uy by the definition of complementarity. 
cP U; = lay by the definition of complementarity. 
4. 


Adding the equations in steps 2 and 3 gives 


Lt by =U; + Us. 


5. Since only even spacings are concerned in the sequence 
quadruple form, when related back to the original (A,B) 
pair, the likes of A are equaltothe unlikes of A for all 
even spacings, This will be expressed as Ly = oan for 
even space, 

6. L,=U, and U,=L, by the definition of complementarity, 
Therefore L,=U, for even spacings, 

ye since the total number of matches at any spacing 

is equal to the likes plus the unlikes, and is dependent upon j, 

10, Adding 8 and 9 gives L 

(U5,,+U 1) = ck; ° 
i. Liyt+Lyy = k; , by step 7, 


jtU +i + ris 


mt ery tU ry = (epytbyy)t 


a 


ee L,,-U, 


in step ll. 


Te O by subtracting the equationin step 8 from that 


13. Therefore Ly, =U nny and the (JII, IV) pair is complementary, 
Since the proof would have been the same if (III], IV) were 
chosen complementary rather that (I, II), the general state- 


ment is true, 


ok oe 3c 
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Theorem 3,2 


If (A,B) are a complementary sequence pair in sequence quadruple 


form and if one of the pair (I,II) or (III,IV) is not complementary 


neither is the other, 


Le 
2. 


Assume (1,II) not complementary. 


There are two possibilities for (IIIl,IV) either complementary 
or not complementary. 


Assume (III,IV) are complementary, then by the Theorem 3,1 
the (I,1IIl) pair must also be complementary, but this is con- 
trary to assumption, 


Therefore (III,IV) cannot be a complementary pair. 


ox os x 


Theorem 3,3 


A necessary and sufficient condition for (I,II) and (UI,IV) to be 


two complementary pairs when a complementary sequence pair is 


written in the sequence quadruple form is if the likes for even spacing of 


one of the code pair is equal to the unlikes at the same spacing for the 


same code of the pair. First the sufficient portion: 

i L, = U, for even spacings, 

ae Lo T Li, = Ly for even spacings. 

on U, +U,, =U, = Ly for even spacings, 

ae L, +L, + U,t U,, = ol, for even spacings by adding the 
equationin step 2 and step 3, 

> Ly 5 U, = ss since the total of likes and unlikes must equal the 
possibilities, 

6. epg eI 

7 L, + U,+4,+ Ur = as - clay by adding the equations in steps 
5 and 6 and comparing to those in step 4. 

8, L BE Li = a rewriting step 2 with ig substituted for L,- 

oe Ly, - U, = 0 subtracting the equation in step 5 from that in 
step 8. 


a 





10, 


es 


ee 


vee 


14, 


i. 


Therefore L,, = U, and (I1,II) are a complementary pair and by 


Theorem 3,1 (III,IV) are also a complementary pair, 


For the proof of the necessary condition assume that (I,II) are 


complementary and that Tis < Li, for some even spacing. 


Therefore U,+U,,< b+ L,, for some even spacing. 


However, by assumption (I, II) are complementary and 


L, = Us, and UL = Lie 


Adding the two equations in step 13 gives Lt+Liy = U;t+Uy, : 


This is incompatible with step 12, therefore it is necessary that 


Uy =L,y or Up = Lp for all even spacings in order to have 
(I, II) and (III,IV) complementary. 

* « kK 
Theorem 3.4 


If the sequence quadruple form of a complementary pair can be 


transformed into the pattern (I II I II), either the original pair is the 


quad or (I,II) and (II,IV) are complementary pairs. 


Bs 


Given that (I II III IV) is transformable into (I Il I Il), and 


is so transformed, 


Lt Lyythyyyt hyve eye yt YU yt UF ag by the supplementary 


property, Theorem 6.1. 


L -br and U, Fur since 2 ave = ae by Theorem 2.1. 
Therefore eL,t+¢Li= 2U/+2U 51> aS by substituting into step 2, 


and using the condition of step |. 


L thie nS : 


L,+U,= Se since the total of likes and unlikes must equal all 


possibilities, 

L,-UF 0 subtracting the equation in step 6 from that in step 5, 
Therefore L,,= U, and (I,II) are complementary. 

Since (I,II) are complementary, soare (I, lI) and any general 


transformations of (I,!I) and (III, IV). 


28 >I oS 
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This last property should be extremely useful in decomposing long 
codes, since it avoids the tedious task of checking to see if (I,JI) are 
a complementary pair, Whenever a complementary pair expressed in 
sequence quadruple form has the pattern (I 11 I TI), certain additional 
transformations are allowed which are not true in general. These 
allowable operations are the time reversal of I and II], II and IV, 
or both, Also allowed as an operation is the altering of either even bits 
or odd bits for all four sequences, The proof of all these operations 
being allowable is quite simple, but does require an addition to the no- 
tation to cover these particular operations, Double or quadruple sub- 


scripts will be used to denote these operations, T 13: an and T 


24 1234 
represent respectively the time reversal of (I, III), (II,IV) and the 
time reversal of both pairs. Als represents the altering of the odd 
bits in all four sequences while A>4 stands for the altering of even bits 
in all four sequences, Actually complementing any two, or all four of 
the four sequences will also leave a complementary pair when reassem- 
bled, but each of these operations is already accomplished by the 
general transformations, 

Mincorerm 3,5 

Whenever a complementary sequence pair is put into sequence 
guadruple and then transformed so that the resultant group pattern is 
felt I Uj), | and Til, or II and IV, or both, can be time reversed 
and the sequence pairs formed on reassembly from thequadruple, after the 
Operations (T) 3 T54 or T1234) will still be complementary sequence 
pairs, 


1. Given, (A,B) is a complementary pair of length n transform- 


able in quadruple form to (I Il I Ii), and is so transformed. 


2. (I, II) is a complementary sequence pair by Theorem 3,4, 


4] 





3. Therefore (I,II), (I,II) and (I,II) are complementary pairs 


by Theorem 2, 3, 


4. From this it is seen that the original sequence quadruple after 


13? +24 and T1234 ane oe lg Ty), 


(I I I Ty) and (I I I Ti). When interlaced in the customary 


the operations T 


manner to form pairs, the pairs they form will be complementary 


sequence pairs by Theorem 2, 8. 
NS oS * 
Theorem 3,6 
Whenever a complementary sequence pair is put into sequence quad- 
ruple and then transformed so that the resultant group pattern is (I II 
ij ay the operation Ay 4? the altering of odd bits, or the operation 
Ad 4? the altering of even bits, will allow the quadruple upon reassembly 


to still be a complementary pair. 


l. By hypothesis A and B are a complementary pair which 
when expressed in sequence quadruple form are transformable 


ter Jiel Il); andis so transtormed, 
2. (I,II) 1s a complementary sequence pair by Theorem 3.4. 
3. (I, II) when operated upon by either Ay or A, is stilla 


complementary pair by Theorem 2,4. 


4, If the symbol (1,2) is used for the pair (I,II) altered itis 
seen that the pair formed from the interlace of (1 212) is 


a complementary pair by Theorem 2, 8, 
x ** S 
These last two theorems should prove quite useful when trying to 
estimate an upper limit on the possible number of composite comple- 
mentary sequences of any particular length. An example to illustrate 
the use of Theorem 3.5 might prove helpful, Taking pair number one 


from Appendix II, 
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A=1101000111011110 
B=1000010010001011 
I=10001011=I1 

Bel Ptoll=1) 

I1J=10001011=1 

IV =1 0000100 =IT_ 


— 


Since this is already in (I II I II) form without transforming it, 


T 3 can be applied directly giving 


I=11010001 

II=01111011 

II=11010001 
IV=10000100, 


Reassembling gives 


Pe) VYOO1101 010110 
B+=1010011000000011 


which is identical with number three in Appendix II, 
As an example, to illustrate the use of Theorem 3.6, applying Ao, 
to the same (A,B) pair gives 


1,=11011110 
2,=00101110 
3,=11011110 
Za o1000] ; 


reassembling give 


A*=1011011110111000 
B*=1110001011101101., 


This is the same as number 6 in Appendix II with the transformation 


w@--© applied. 
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Unorde red Unordered 


Transform Equivalent Transform Equivalent 

l. i= 1 1 Ii) ies jee es dil age O 
ae fea) il I iy C 18oaO =e 1 a J 

3, Peal al I 11} I 19, Pez=(III I J) si 

Ae ie elel 1 IT) S 20,2) Ov aee ata R 
5. ie iee Il I) xX ei)6CUWNN = (1 ee 1) U 
6, ete ir 1) P 22.5 ENE (a Seis 6 
Te @e=(1 .i1 1 If) F 23, 0 = 1 ie V 
8, Ce=(I II I I) D 24. Kz=Q0T1 I 1) Vd 

9. €e=(1 II I 11) a 25. =) eels Q 
10, Z =(II1 T T TM)) K 26. H=(Il lu y W 
ale Y=(1l I 1 Wy) M 27,  G = (ieee tien) B 
lee, X =(11 II It) To 28. Fe=(I I I Ip Cy 
13, W =(II I I I) H 29 D = (Ue Her oil) Ge 
lea V =(I II It Jj iu 30. we(ll In yj .) 
15, U =(I II II 1) N 31, Bei iy G 
16, S =(I II II J) Ten 32500 O = (ile ieee) TT 

= (ie ieee Ao, I 


Redundance due to (I II I IJ). 


TABLE 3.4 
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fire stioIt I TI) form, or one of its transforms, appears quite 
often in the study of composite complementary sequences, Since two 
of the components are the same, and the other two are complements, 
it is reasonable to assume that some of the transforms are redundant, 
that is yield indistinct results. Table 3,4 lists all the transforms 
with the elements I, II] and their complements, If (A,B) are taken 
as an ordered pair, the number of different codes formed from the 
transformation group are 32 instead of 64 because Agzg(I II I TI) 
and E(I II I Il) are the same. If (A,B) are taken as an unordered 
pair the number of transforms forming different codes is decreased 
from 32 to 16. Table 3.4 lists each transform with its unordered 
a be 

The only possibility of eee less than 16 unordered pairs of 
codes which are distinct would be if the group form could be transformed 
into the pattern (I I I I), The only case known to the writer where 
this is true is the quad. For one version of the quad this is 


A=11 B=10 
l=] 
II=] 
IlI=1 
IV=0, 
Table 3.5 shows the eight ordered pairs, or if the cross index column 


is used the four unordered pairs, for the quad. The four possible code 


transformations for the quad are: 


1 vs 3 4 
A= 11 Lied 00 00 
B= 10 01 10 O01. 
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Transforms with pattern Cross index number 


Pattern for ordered pairs for unordered pairs 
eet 1 Td) I,T1,H,t 3. 
eel 1 Y 1) Ai, Y,V,B 6. 
emer (i I 1 1) Az, W,S,0, E l. 
te {t i 1 iO) fe EOE IN 8. 
eed I Ij CO, a 
oat 1 I Xj Ce N2,.G a 
feeee(t I 1 Jj Ties nO A a 
ome) I I J) X,U,P,D 4, 


Redundance due to (I I I J), 
TAB EE AS. 


This exhausts the transformations of special types and the redun- 
dancies caused by (I II I II). The operations group will now be 
compared to Senior's group number 44 of length 32 as was mentioned 
earlier inthis chapter. 

The identification of any group with a group in Senior's list, is 
accomplished through establishing a one to one correspondence between 
eener ators, As will be demonstrated, Senior's group num- 
ber 44 of order 32 is in one to one correspondence with the unordered 
Operations group of complementary sequences expressed in sequence 
quadruple form, Therefore the two groups are isomorphic. 

Senior's generators will be given in lower case letters, while those 
of the operations group will be in their standard upper case symbols, 
Five generators are required for Senior's number 44; let these be 
a,b,c,d,e. The following relationships between generators are re- 
quired for Senior's 44: 


i? a® =b2=c%=d*=e" =1 where 1 is the identity. 


at 


ae b Ab=a, c ltac=za, d ‘adz=a where x ~is the inverse of x, 
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Be ey De=p d. tad=b, 
4. d 4cd=c, 

oe eu ae=a, e ‘be=b, 
6, e lce=ca, e 'de=bd. 


A check of the above requirements shows that a,b,c,d form an 
abelian group of order 16 while the "e'' generator is commutative 
only with a,b andis non-commutative with c,d. 

Let a=C,, b=Csz, c=T,, d=Tz, e=A,,. From Table 3.3 it is seen 
that the following relations are true, 

eee = C5= 17=Ts= Af=l 

ae Gre, Ga= Ci, -1T7 Ci TiS Ci) 12 G =] Ce 

a: Weel, - C2. 

4, Wells = 91). 

oe AiG; Ai, = Gi, Ai Cp Ai =Cos., 

6, Ai T, A, = Ti Ci,= 0, A, Tz A, = Ca Tse H, 

This one to one correspondence between the generators of Senior's 
group number 44 anda set of generators from the unordered operations 
group on complementary sequences is necessary and sufficient to prove 
the two groups are isomorphic, 

The diagrams for Senior's group number 44 of order 32 and 
those for its subgroups were used as a guide line to check the generation 
of subgroups from the operations group of unordered complementary 
sequence pairs. These subgroups of order 1,2,4,8 and 16 are shown 
in Appendix IV, All groups that were shown in Senior's diagrams are 
included in these charts; however, Senior's diagrams included only 
normal subgroups and, although much care was taken in the generations 
of the subgroups of the operations group, there is a possibility that 


Appendix IV is not exhaustive due to some undetected non-normal sub- 
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Senior in conjunction with others has rewritten his work and 
changed his designating method. ne Under this new method his 
number 44 of length 32 is now May. The ordered operations group 
of length 64 is also identified as [55476 

An identification of the group with Miller's number 23 of order 
32 was also accomplished for those who prefer the substitution group 


method of identification, Le 
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CHAPTER IV 
HAMMING DISTANCES OF COMPLEMENTARY SEQUENCES 


The Hamming distance is a well known property of certain classes 
of codes, a This very useful property has had most of its applica- 
tions in the fyelds of binary error detecting and error correcting codes, 
This chapter will demonstrate the usefulness of Hamming distances in 
the field of complementary sequences, as both the code pairs and their 
sequence quadruple formulation have certain invariances in Hamming 
distance under various conditions, 

The Hamming distance, D(U,V), of two binary vectors or sequences 


U and V is defined as the number of positions in which these two binary 


vectors differ. 


1=n 


D(U,V) = u, @ v, (4, 1) 
i=l 


For example, the Hamming distance of the two binary vectors 


U=l1011001 
V=0111001 


is D(U,V) = “ u; © vee (1+1+0+0+0+0+0) = 2 
i=] 
The Hamming weight or just weight of a binary vector is its distance 
from the null vector. 
Theorem 4.1 
All complementary sequence pairs of length n have a Hamming 


distance of n/2. 


By hypothesis 


A and B forma complementary sequence pair, 
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le, All complementary sequence pairs satisfy the parity check 


a. ait] ®b. 2) oR | l, for all 
eles i <= 1/2 by equation 2, 4, 
le Since mod 2 addition is commutative 
ae eee ae Cnet ame 
iF Therefore either (a, ®b.) = 1 vor (a is4y @b _s41)? but not 
both, 
4. Foreach i, 1¢<ign/2 one of the two vector positions in 


step 3 will have a one and the other vector position a zero, 


Summing up overall possible i, gives n/2 for the 


Hamming distance. 


When a complementary sequence pair is decomposed into a sequence 
quadruple, each pair of sequences in the quadruple, which form one of 
the complementary pairs, has a Hamming distance which will be shown 
to be characteristic of the code length. The sum of the two pairs of 
Hamming distances from the quadruple is an invariant and is equal to 
n/2 the Hamming distance of the complementary sequence pair. 

For example, 


A=1001010001 
B=1000000110 


are a complementary sequence pair of length 10. 


Decomposing them into the standard sequence quadruple, (I II III IV), 


gives 
J = 10000 
Ii = 10110 
III = 10001 
IV = 01000 
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I 9 II = 00110 D(I, ll) = 2 
III 9 IV = 11001 D(III, IV) = 3 
I @ Il @ Il] ® IV=11111 D(P@II, III@IV) =5. 


This distance is half the length of the complementary sequence pair, 
and is required because of the necessary condition of the parity check, 
If the distances (2,3) forthe pairs given in the above example are 
taken as an unordered rather than an ordered pair, this (2,3) pair 
will be invariant through all general transformation operations, 
All of the general transformations on the sequence quadruples 
have one or more of the following properties: 


le Change the order of (I,II) or (III,IV) or both, 


for example, T, , Ta, iia 


oe Complement the pairs (J,II) or III, IV) or both, 


for example, C, ,C2,C. 


on Complement one from the pair (I,II) the other from 


the pair (IIJ,IV), for example, A, , Ag. 

4, Exchange (I,II) for (III,IV), for example, E. 

Since all operations in the group can be generated by multiplication 
from the set listed above, an examination of the effect on Hamming 
distance by these transformations will be adequate to prove the invar- 
lance, 

Theorem 4.2 

A change in order of (I,JJ) or (III, IV) or both will not change 
the Hamming distances of the pairs (I,II]) or (10, IV). 


n 
lV By definition D(U,V)= 2 u. ® v. 


1 
i=l 


Le u. @v. = v. Gu. modulo 2 addition 1s commutative 
i a 1 
ral rea 
ce Die == | > u, Ov. = L2 V; Gu. = Dey; un 
ak I T=] 
% 26 *K 





Theorem 4,3 


Complementing (I,II) or (III,IV) or both will not change the 


Hamming distances of the pairs (I,II) or (III, IV). 


he Tow. = G. & V. by Theorem 2,1 
1 1 1 1 
ak Pherefore yu. vy. = = vu) 687. 
1 1 1 1 
ce and D(U,V) = D(U,V) 
K xk 6 


If it is assumed that the length of the sequence quadruples formed 
from the complementary sequence pairs have a length K and that the 
distance of the first pairis R, then the following theorem can-be proved: 

Theorem 4,4 

Complementing one of whe pair (I,1l) and one of the pair (iia 
results in the exchange of their Hamming distances, or expressed dif- 
ferently, gives the same distance pair in opposite order. 


Le Let D(I, Il) = R, therefore 
D(III, IV) = K-R 
as the parity check requires D(I@II, III'@IV)=K 


Z, Complement one of the first pair inthe (III II IV) form, 
and also complement one of the second pair, 
Each vector position that was formerly alike and summed 
to zero is now different and sums to one, 
Similarly each vector position that was formerly different 
and summed to one is now alike and sums to zero, 


BM D(I, 1) = D(I,ID = K-R from 2. 
4, Dili IV) = D(ill,IV) = K-(K-R)) = Betrem 2, 
* 1e * 
Theorem 4.5 
The exchange operation has the effect of changing the order of 
Hamming distances of the pairs (I,II) and (III,IV) when expressed 


in sequence quadruple form. 
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Fe The pairs (1,II) and (IIJ,IV) are not changed in any way 
except by position in the exchange operation, 


rap The only way to effect the Hamming distances of a pair of 
vectors is to change some of the components of the vectors, 


Sh Therefore the two weights are the same but reversed in 
OFGeYT . 


An examination of the Hamming distances for the known kernels of 
2, 10 and 26 is now in order as these distances are a characteristic 
of the kernel, as are the vectors formed in the process of obtaining the 
distances by the modulo 2 sums, | 
The first case to examine is the trivial case for the kernel of 
length 2, 
A=11 B=10 


Broken down into sequence quadruple representation 


I=] D(I, I1)=0 
II=1 

III=1 D( Ill, IV)=1 
IV=0 


There are two kernels for length 10 to be examined, 


A=1001010001 
B=10000001 10 


I=10000 I @ II=00110 D(I, I= 2 
II=10110 
hia FOooT III @ IV=11001 DC = 
IV=01000 


and 


A=0101000011 
B=00001001 10 


I=00001 1@ II=10010 D(I, II)= 2 
W=10011 
I1I=00101 III ® Iv=01101 D(IiI, IV)= 3 
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Both kernels of length 10 have the unordered Hamming distances of 
(2,3) when expressed in sequence quadruple representation. However, 
their ''Hamming vectors'', the unordered pairs (00110,11001) and 
(10010,01101) are different and are characteristic of the kernel of 
their origin under all possible transformation operations, Hamming 
vectors will be considered at length in the next chapter and proof of the 
above statement will be deferred until then, 

There is only one kernel of length 26 (Chapter 7) and it is 


A=01001101111010111100111010 
B=10110010000111111100111010 


I=0010111110111 I ® II=0000001101010 D(I, I) = 4 
II=0010110011101 


ITI=1101001110111 III ® IV=1111110010101 D(III,IV)= 9 
IV=0010111100010 


This is a complete listing of the Hamming distances for all known 
kernels and it is exhaustive for n=2,10,18 and 26, It is worthwhile 
to note that the sum of the unordered pairs of distances is always equal 
to one half the length of the code and that the Hamming vectors of the 
pairs are complements. Both of these conditions are due to the parity 
check being a necessary condition. 

A natural extension of the Hamming distances of kernels is the 
consideration of Hamming distances of composite codes when put into 
sequence quadruple form. Six different code pairs of length eight are 
listed in Table 4,1, These pairs will be examined for a possible in- 
variance. After this invariance is noted a general proof will be given 
for the characteristic distance of composite codes, 

The first four examples in the table when decomposed into sequence 
quadruple representation have (I,II) and (III,IV) as code pairs which 


are still complementary. These codes are half the length of the original 
pair, and being complementary, have a Hamming distance one half of 
their length, or n/4 as compared to the original codes of length eight, 
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(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


Ae1T111001 
B=11001010 


=PEro 
M=1,0:1 | 
Maio) 
IV=0001 

A=11110110 

B=11000101 


I=1101 
Ii=O111 


fet=1000 

v= ULO1 
A=11101101 
B=11100010 


I=1110 
m= TTor 


ItI=1101 

IV=0001 
A=11011110 
B=11010001 


t=O ael 
=O tT 


IITJ=1000 

IV=1011 
A=11101011 
B=11100100 


T=1111 
Tl=1011 


IITI=1100 

IV=0101 
A=10111110 
B=10110001 


J=l1l11 
II=0110 


IITI=1100 
IV=1010 


H=(00110011) 


I ® IT=0101 
III @ IV=1010 


H=(00110011) 


I ® II=1010 


III ® IV=0101 


H=(00001111) 


Pe I 00 10 


III @ IV = 1100 


H=(00001111) 


I ® II = 1100 


III ® IV = 0011 


H=(00001111) 


I @ IIT=0110 


III ® IV = 1001 


H=(00001111) 


I ® II=1001 


III ® IV=0110 


Dl, 1l)#s"2Z 


D(III, IV) = 2 


D(I, Il) = 2 


D(III, IV) = .2 


D(I, II) = 2 


D{ Illi eae2 


D(I, I) = 2 
lldat ny) = 2 
DUA g—ac 


D(III, IV) = 2 


D(I, II) = 2 


D(III, IV) = 2 


Some Hamming distances and vectors 


of codes of length 8, 


TABLE 4,1 
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An examination of the last two examples in the table shows that 
although the (I,II) and (III,IV) pairs are not complementary their 
Hamming distance is still n/4 when compared to the original codes 
of length eight. That this feature is true in general will now be proven 
by examining and exhausting all the known ways of generating composite 
codes, 

Theorem 4,6 

Every composite code of length n (generated by one of the standard 
methods) when expressed in the sequence quadruple representation will 
Mawee OCT). = D(III,IV) =n/4. 

There are presently known four general ways of gene rating com- 
posite codes from shorter codes, (there is also the method applicable 
only to codes of length 2* given at the end of the proof), These are: 


1. Time Sequence 
5S, = AB 
Sis Saul) 

2. Interlace 
ea 
T2 = 


eo ne 


b,apbs..-a 


n 


ol 


ay 


ele at 


3. Time Sequence Exponential 
Cy Co Cc d d> d 


U,=A A ey ie): 'B eee) 
d d a. =e Cc Cc 
(eee Ar!) Ap eee 
4. Interlace Exponential 
Se nee, m 


Cc 
Te =A GB AB .....A o8 ; 
C dlieegc 


d Cc d 
Vee AB OM m-lyom-l a tp? 


fate 


Throughout this proof it will be considered that (A,B) are a comple- 
mentary pair of length r and (C,D) are a complementary pair of 


Femeth sii. 
S16) 





ihesearity check a. ® a @ b. € b = | is a necessary 


ntl -1 Nt loi 


condition for a sequence pair to be complementary. One of the simplest 
methods for making the parity check is to fold both codes in the middle 
and double them back on themselves, ? This procedure is shown in 


Pigure 4.1. Each 


oA Wear eae, ©ecoceee at a ee 
mena oe et Ontl-1 a a 

rere ny /24+Pn/241 
DID aisle ete b. J22 ope 
Phar yaaa Bat t-i °° Pnf2t2n/2tl 


Foldover Parity Check Method 
Figure 4/1 


column in Figure 4.1 satisfies the parity check. If code pairs are 
written in the sequence quadruple form suchas in Figure 4,2, itis 
seen that the columns of this array are the same as the column in the 
first array with rearrangements in order provided it is allowed that 


the columns 


a a 
ic u 
a and a 
u r 
b b 
r u 
b b 
u r 


are the same. Since modulo 2 addition is the Operation used in 

making the parity check and in determining Hamming distance, and the 
use of each column involves only modulo 2 addition which is commu- 
tative, the two columns can be considered the same. Therefore the 
fold-over method of grouping and the sequence quadruple form of grouping 


have invariant columns, order being ignored. 
oa 





auc sae aa 
ec cee a, ap 
BERL te eee: bab 
DADE. >. aan 


Sequence Quadruple Form 
of Parity Check 


FIGURE 4, 2 
Now consider the composite generating methods, 
1. Time Sequence 


$4 = AB 
S2=AB ~~ writing in the fold-over parity check form 


ee Do 3 eas nae a n= 2 fee 


D(A,B) = r/2 
> Sy cobs Ba Wn 8 Dy 


oo aa BYU Jsiuee ee 2 
r r-l 3) Z l 

The top pair are complementary and therefore have a distance of 
r/2, similarly the bottom two rows are complementary and have a dis- 
tance of r/2 by Theorem 4.1. To write the time sequence code in 
sequence quadruple form would not change the distance since it would be 
purely a reordering of the parity check columns, Therefore D(I, II) 
would equal r/2 as would D(IU,IV). Interms of the composite code 
of length n D(I,II)=D(III,IV)= r/2 = n/4 whether or not (I,II) and 
(III, IV) happen to form complementary pairs, 

2. For the interlace scheme the (I,IJ) pairand the (II,IV) pair 
are obviously complementary sequences by construction and will have a 
@ictance of Ly 2 = n/4, 


3. Writing the time sequence exponential pair in fold-over form 


eats € 
gives ween Ly 4 am 
q where the symbol B means B time 
B ae aS ewee B 1 reversed, 7 
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The B's are all time reversed but the time reversal of one of the 
complementary pairs does not change the complementary property 
D(A,B)=D(A,B) = r/2,. For each bit of the C code and D code there 
is a distance r/2 since every A is matched with some transformation 
of B for each of these bits, There are a total of m bits inthe C 
code and inthe D code, therefore the total distance of U, folded back 
on itself is mxr/2=mr/2. Similar reasoning gives U2 folded back 
on itself a distance of mr/2. A rearrangement of the columns will not 
change the Hamming distances Therefore D(I,ijJ=mr/2 and 
D(IlI, 1V)=mr/2, since the total length of the code is n=2mr, D(I, II) = 
nil, 1V)=n/ 4. 

The interlace exponential scheme is the last general method of 


code generation to be considered, Written in the fold-over form it gives 


Cc d c d e d 
en le OK. 25 2 ga mlep m/2 
d e d d c 
B av FS m-1 _B Ca m/2+1 
d C d d Cc 
Wee Mp My, mo, A m/2+1, m/2+1 
Cc d Cc Cc d 
B la cai B m/2 Agen 


Again each A is matched with some transform of B, D(A, B)=D(A, B)= 
D(A,B), etc. = r/2. As before, rearranging columns into sequence 
quadruple form will not effect the Hamming distances, The distance of 
the top two rows is mr/2 andthe bottom two rows (Vz) also is equal 
to mr/2. The total length of the code is 2@mr, therefore 
Die jill, 1V j=mr/2=n/4. 

oe When composite codes are of length hae they are sometimes 


formed ina special way from the interlaced sections of a complementarv 


oh, 





pair. These sections or pieces of the codes must be of length p=2) 


ee Z = 
where j<r. Using the symbols an a) 4a,aza4--- -2.53 i Pe ratio 4 


220 and so, the following pair would be complementary if A and B 


formed a complementary pair. 


pre eA Bo. oe ts 


Is B 
lp pp 2p 2p np/p’ np/p 


Sa SN ei Ney aN B 
2p p p 2p 2p np/p np/p 

Rather than try to use the general term to prove that D(I,II)=D(III, IV) 
=n/4 as was done in the previous cases, a short general example of 
length 8 will be used. These two codes (A,B) of length eight will be 
interlaced four bits ata time (p=2) to form a code of length 16, 

Mes 12233 a 2 8452627789567 8 

S 54742723446 )b,b4b,a5apa7agb5bgb 2b, 
Writing this composite pair in fold-over form yields 

Oa eh ae op Don 

bob-b b,a,arayae 
Spa 2 Goan pees 


Brees ce7 7625 
D(A,B)=4 by Theorem 4, 1. Sia folded over is just a column 
rearrangement of the complementary pair (A,B), and sincea re- 
arrangement of columns does not effect Hamming distance, the Hamming 
distance of Sig folded-over is also 4. Similar reasoning holds for S04. 
Since the standard sequence quadruple form is again just a re- 
ordering of the columns, the distance of the (I, II) pair and the 
(III,IV) pairis n/4, For this example, 
I=a azb,b,a,ajb,b. 


: D(I, I)=4 


Il=bebpagapb,b,a4a> 
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Ill=a,a.b,b b 


b 
hes 3 5° 75 
D(III, IV)=4 


ibe b,a,a,b,D>5a,a 
2 
TE a uC MT 


Therefore any composite complementary sequence pair when written in 


sequence quadruple form has D(I,II)=D(III, 1V)=n/4, 
ae aK xe 


Although there is no proof that the 4 standard generating methods, 
plus the special one applicable only to ee codes, are the only way to 
form complementary codes which are composite, an exhaustive search 
of code lengths 16 and 20 revealed no codes which were not formed 
by these standard generating methods, 

The full usefulness of the’ Hamming distance property is, of course, 
not yet known, The property of the distance of the complementary pair 
being equal to half the code length as shown in Theorem 4.1, is used 
in many of the proofs in the chapter on Hamming vectors, This property, 
D(A,B)=n/2, indicates that the code pairs are orthogonal, Because of 
this orthogonality property the leakage of an A code carrier into the 
B code filter, or vice versa, will have no effect on the output at exact 
match, 

The primary use the writer has found for the Hamming distance of 
the (I, II) and (III, IV) pairs, has been the error check that this invariant 
property offered in decomposing codes by hand, One possible use for the 
invariance of D(I,II) and D(III,IV) could be in the search for kernels of 


new lengths. If a characteristic Hamming distance for unknown kernels 
could in some way be deduced, based only on length, this invariant prop- 
erty of distance would act as a most powerful screen in the search, The 
(I, II) and (IIl,IV) pairs are themselves orthogonal, unless the code is a 
kernel, This may provide other possible methods for searching for new 
codes. It is difficult to say in advance what utilization will be made of any 


invariant property, but it is fairly safe to say that the property will find use. 
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CHA PIER 
THE HAMMING VECTORS OF COMPLEMENTARY SEQUENCES 


The Hamming distance between the code pairs in a complementary 
sequence was shown by Theorem 4.1 to be one half the length of the 
code. This Hamming distance was obtained by first finding what I shall 
call the Hamming vector, H(A,B)= A@B or for convenience just 
H=A®B, and then counting the number of ones H contains, The Ham- 
ming vector is itself useful in recognizing the composition of a composite 
code when it 1s desired to break the code down into shorter lengths. This 
decomposition might be used to find previously unknown kernels, or to 
set an upper bound onthe possible number of pairs of codes of any length. 

An examination of each of the methods for generating longer codes 
from shorter ones as explained in Chapter 2 will show that each method 
has a characteristic Hamming vector form, although it may be necessary 
to time reverse one of the codes to obtain this form, The method of 
generating codes of length 2, which does not apply to complementary 
codes in general, will also be studied with its characteristic Hamming 
VEGtOg, 

In the following proofs (A,B) form a complementary code pair of 
length n while (C,D) form a complementary code pair of length m. 
The symbol (0)” or just 0” will mean a string of n zeros, similarly 
(01)" will symbolize a string of n zero ones, and (a.@b.)" will be a 
string of n elements equal to (a. ®b.) : 

Theorem 5,1 

The time sequence form of generating a composite code has the 
characteristic Hamming vector H=(0)"(1)", 


Use the time sequence generating method on (A,B) to form 
5 = AB 
S2= AB . 62 





Take the Hamming vector of (S,, Se), H=S; 6 S2= Nee this is 
(50) 


seen to be true by inspection and definition of modulo 2 sum. 


Xe 2 oe 
Theorem 5.2 
The interlace form of generating a composite code has the charac- 
teristic Hamming vector H=(01)”. 
Use the interlace generating method on (A,B) to form 
T\=a)b,a,b,a,b,. AGAR AE abo 
Ta=a,b,a,b,a,b PRR: 5c ekeUI e 
Take the Hamming vector of (Ti, T2), H=T, T= Ollie (3.2) 


this is seen to be true by inspection and definition of modulo 2 sum. 


oe oe *K 
Theorem 5,3 


The time sequence exponential form of generating a composite code 


has the characteristic Hamming vector H=(c,@d ye (c,8d ye 
er ara 2 mM 


-1 
n = n = n 
(c_@d,) (d, @c_) ; -..d_ Oc) ) : 
Use the time sequence exponential method on (A,B) and (C,D) 


to form d 4 


re @ Cc d 
U,=A A ee ete eer: Aes SF oA 
Cc 


B 
dd e c 
Us=A SI. m-1 e@e#e8¢e@ 8&¢ 6 86 @ re B gras m-1 9 8 @ @@¢ &© @ B ; ® 
Take the Hamming vector of (Ui , U2), 


n n n = 
H=U,@ Us =(c,@d_ ) (c,8d__ ) voeee(c Od)) (d, Gc ) 


-| 
ears ( 559 
(d_6 Cc) : 


This is seen to be true by inspection and definition of modulo 2 sum. 
ot 3c a 
It is to be noted that each of theTheorems 5,1,5.2,5,3 have 
Hamming vectors which are antis:ymmeric about their center, This is 


required to satisfy the fold-over method of making the parity check, 
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Theorem 5.3 also shows that along with the anti-symmetric characteristic 
the Hamming vector for time sequence exponential composition has unit 
clusters of n zeros or n ones and is also anti-symmencin both halves, 
This is easily seen because the first sum is (c $d _) while the last sum 
to operate onthe A code is (cd, ), but it is known through the parity 
check that c,8d_ Oc _ Od, = 1. Therefore one of the sums (c,4d_) or 
(cod, ) must be zero and the other must be one, To summarize, for 
the time sequence exponential form of composition the first half of the 
Hamming vector is antisymmetric about its center (the quarter length 
point) and the first half and the last half of the vector are also anti-sym-+ 
metric, 

Theorem 5.4 will concern the interlace exponential method of 
forming composite codes and it will be seen that this method will again 
form a Hamming vector anti-symmeticabout its center, but with a dif- 
ferent arrangement of the clusters of zeros and ones, 

Theorem 5,4 

The interlace exponential method of generating a composite code has 


— | 7 n al n 
the characteristic Hamming vector H=(c)4d_) (dc) (c,@d__)) ere 


n =e ol 
(c_@d,) (d_8c 1) 
Use the interlace exponential method on (A,B) and (C,D) to form 
Cy d, C5 d, cm dn 
Vi=A lo peat TP eect bie A B 
~m .~m-=-1 oF Sul 


Take the Hamming vector of (Vi, Va) 


n - ,n n n n 
H=Vi, 9 V2 (c)8d_) (d) 8c) (c58d___)) oeee(c Od, ) (d_8c,) (5. 4) 


This is seen to be true by inspection and definition of modulo 2 sum, 
2c xe * 
The interlace exponential method again gives a Hamming vector 


which is anti-symmetric about its center due to the condition of satisfying 
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the fold-over form of the parity check, but its clusters of zeros or ones 
are now én in length rather than just n as was the case for the time 
sequence exponential method of generation, This characteristic of 
clusters of 2n in length is caused by the necessary condition of the 

parity check for the (C,D) code pair, For example, the first two clusters 
are (c,9a_)” and (d,@c_)” and the second two clusters are (c,@d_ 4)" 
and fabio)”. Applying the parity check to the i=l and i=2 bits 
in the (C,D) pair gives the following two equations: 


© Oc 0d od _=1 ey ante c5Oc__ 9d,0d_ = ; 


] 
Adding one modulo 2 to both sides of both equations and rearranging 
gives (c )$d_) O(d 9c 1) =0 and (c,@d___)) O(d,@c__ 1% 1)=0, Itis 
true in general that x ® i = x; therefore (c,d _) ®(d,Oc__)=0 and 
(c,@d__)) ® (dj _,) =0. The sums (cy @ dn) and (d) ® cay) must 
be alike to satisfy this constraint with each of the sums giving n ones 
if unlike, or n zeros if alike, and since the pairs are in time sequence 
the clusters must be in groups of é2n, A similar reasoning will show 
that the parity check for i=2 gives identical results for the third and 
fourth terms in the Hamming vector for the interlaced exponential form 
of generating composite codes. Therefore equation 5.4 can be written 


2n 2n - ,2n 


H(Vi,V2)=(c,@d_) > (c, @d_ 1)... (d_ Oc)) (5, 5) 


Before going into the special case applicable only to codes of length 
2” it might prove helpful to give an example of generating a code by one 
method and then decomposing it, to show that it could have been composed 
by a second method. The time sequences exponential scheme will be the 


generating method used, 


A=10 C=1001010001 
B=11 D= 10000001 10 
U,=1001011001 10010101101 1000000000000111100 
Us=011010010101010101 1000111111001100111100 


H=U,®Uez =1111111100110000000011111111001100000000 
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In passing it should be noted that the shortest clusters of zeros and ones 
in H are of length 2 and that H is anti-symmetric about its center. 


If U, is now time reversed, indicated by UL , the complementary 


pair 1s now: 


U,=0011110000000000001101 10101001 1001101001 


Uzg=0110100101010101011000111111001100111100 
H = U) ® Uz=0101010101010101010101010101010101010101 = (01)*" 
This emi of the Hamming vector indicates that H is composed of an 
interlace pair of length 20. Decomposing both UL and Us by assum- 


ing an interlace pair gives 


U , =01100000010111010110 
01100000011000101001 

H = 000000000011111f1111 = 0/"1'9 

Us= 01100000010111010110 
10011111100111010110 

Mem 1oceodoeeen = 1'°o,” 


The Hamming vector of U, indicates that this code of length 20 


] 
was composed by the time sequence method from codes of length 10, 


If the form AB is assumed for the codes of length 10 then 


A=0110000001 
B=0111010110 


This (A,B) pair is seen to be the transformation ETC,= EG onthe 
Original kernel of length 10 which was used to generate the code bya 
different method, Similarly the pair from Uszis also a transformation 
of the original (C,D) pair. Therefore the given code of length 40 
could have been generated by two different methods, 

A second example will now be given using the same (A,B) and 
(C,D) kernels as were used inthe first example, but his time utilizing 
the exponential interlace scheme of generation rather than the time se- 


quence exponential form. 
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A=10 C=1001010001 

B=11 D= 1000000110 
V,=1011010001001000010010000100011101111000 
Ve=CLOCIOILIOIIOIIIOLOOOL TOL C rr ror ooo 

H=Vi@V2q  =1111111111111111000011110000000000000000 


Note that the smallest cluster is four, as would be expected of an inter- 
lace exponential scheme with the (A,B) pair of lengthtwo. The H 
vector is again anti-symmetric about its center, Time reversing Vj, 
gives 


V, =0001111011100010000100100001001000101101 
V2 =0100101110110111010001110100011101111000 


H=V | ,@Ve =0101010101010101010101010101010101010101 = (oles 
As before this form of H indicates an interlace code pair, Decomposing 


both vi and “5 into interlace code pairs gives 


V , =00111101000100010110 
01101000010001000011 

H =01010101010101010101 = (01) 

Va= 00111101000100010110 
10110111101110111100 

H = 10101010101010101010 = (10) 


10 


10 


a and Vz both break down into interlaced pairs of length 10, The top 


row of FT breaks into A=0110000001, B=0111010110 which is the 
moitetormation EG operating on the original (C,D) pair as in the 
previous example. Each of the codes for V2 would break down in 
similar fashion to some transformation of the original (C,D) code pair. 
One use for the method just shown could be to search for larger 
kernels than those currently known, An example might be kernels of 
length fifty found through the use of time sequence exponential or inter- 
lace exponential, where both (A,B) and (C,D) are codes of length 10, 


These would form codes of length 200. From this as a starting point, 


a decomposition to codes of length 100 and thento length 50 might 
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be possible through this technique, This technique will be discussed 
much more thoroughly towards the end of this chapter, but it is neces- 
sary to consider the special case of 2" code composition and the 
Hamming vectors of the sequence quadruple form before continuing with 
this interesting possibility. 

The special generating method which is applicable only to codes of 
length 2" is actually a generalization of both the time sequence and the 
interlace methods. As mentioned inthe previous chapter, this method 
takes a code pair of length 2" and from them forms a code pair of 
length De . The allowable combinations are formed by interlacing 
pieces of the A code andthe B code of length 2™™, where f=, 0; rae eae 

As an example a code of,length 8 can be used to generate codes of 
length 16. Inthis case r=3 and m has the possible values 0,1,2, 3. 
These values of m give section lengths of 1,2,4,8 respectively. 

A=a)a54,a,a,apazae 
B=b)b,b,b,b,-b, bob, 
fOr Sim=0', 


S 1 72 Pa 2b 24 3b3a4b4acboapbpa lb 2a gb, 


12-26243b,a,byacbpa bya 7b7a gb 
8 1 8 


1 
= 0101010101010101 =(01) =(0 1 ) 


Se =a b 


Zi 


H=S) 95) 


This is exactly the same as the interlace procedure which is a standard 
generating method, 
The section lengths are 2 bits long for m=l., 


S 27a 1a2b)b,a,a4b,byapagb.bpazagb lb, 
S59=a ;a2b)b,a,44b,b,a, 7gP7? 2 
= 0011001100110011 = (0°1°)*, 


a b, bea a (hb 


H=S) 79555 
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m=2 gives 


aya bob. b-bDeaee bebe 


S 34 4h 1 92535,454,4 7a gd bp bob, 


[aa 1-2 


2474 1 4723a,b)b,b,b,a,a,atagb, bp babe 
AeA. 


H=S,,95,, = 2000111100001111 = (0°17), 


S 


For m=3 the lengths are eight, which is the entire code. 


5) g=AB 


S59=AB, 


This is the same as the time sequence form giving a Hamming vector 
H=0000000011111111 = (0°1°), 

The Hamming vectors from this special construction are altering 
clusters of zeros and ones, each cluster being the length of ene, segment 
used in the generation, 

Theorem 5.5 

Complementary sequences of length pea formed from an (A,B) 


pair of length n=2° using the special generating method with segments 


2n/m 
of length m have a characteristic Hamming vector H=(07"1™) 


Use the special generating method utilizing segments of length m 


of the pair (A,B) to form 


b 


Pt 223° mn P1P2°* > Om Amtl* mt2.,.°2m m+] ate 


So ne 2232 8 em Pp Paee OL Aang mde? ' 22m Pmt Pmst2** On: 


Take the Hamming vector of (Syn? Sou 


eam rn, Canaan 
H=S,_,9S,, = (0 1) , 


This is seen to be true by inspection and definition of modulo 2 addition, 
a 4 a6 
It now seems suitable for completeness to include a theorem which 
has been alluded to many times but not proven, This is the anti-sym- 


metric property of all Hamming vectors of complementary sequences, 
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Theorem 5.6 
The Hamming vector of a complementary sequence pair has a pat- 
tern which 1s anti-symmetric about the cénter of the Hamming vector. 


ie A=a)azaza, Rei eee a a. 


n-lon 
B=b,b,b3b, Sere ceeencctre bet |? is a complementary pair of 
length n 
2. H(A,B)= h hoh.h,. = he /2: a ho by definition of Hamming vector 
ole a. @a_,,_6 b. 9 bi = | by equation 2,4 
4. a @ b. = h. 
and both by definition of h. 
Berane: © Pnti-i > “ned 
ae h. @ hot)ei =] 
7. Add (hoa i ® 1) to both sides of the equation 
eaien 7 i 


Make a change of reference to the center of the vector, which lies 
between m/2 and m/2+1, rather than the end by having i=n/2 -r then 
ntl-i=n/2t+l +r 


=h 


Ale 4l+r for all 0 <r<n/2 -1 by step 8, 


oe De) = 
10. Therefore the Hamming vector is anti-symmetric about its 


center when it is formed from a complementary pair, 
* oe ae 

The Hamming vectors of (I,II}) and (III,IV) in the sequence quad- 
ruple form are invariant under transformation if H(I,II) and H{(III, IV) 
are taken as an unordered pair. The following four theorems which 
parallel Theorems 4.2,4.3,4.4,4.5 on Hamming distances will dem- 
onstrate this characteristic of unordered Hamming vector pairs. 

Theorem 5,7 

A change in the order of (1,1]) or (JIJ,IV) or both will not change 


H(I, II) or H(III, IV). “6 





i>. nu, V)=u, @ Vv; by definition, 


Die u; ® een ee ® u, since modulo 2 addition is commutative, 


a. theretore H(U,V) = H(V,U) 
Xe oe x 
Theorem 5,8 
Complementing the pair (I, II) or (III,IV) or both will not change 
Piet) or H(III,IV). 


i mtu, V) = u, 9 v. by definition, 
a geome = UW. Oy. by Theorem 2,1, 
i 1 i i 


oe therefore H(U,V) = H(U,V). 
oe * Xe 
In the following proof it is assumed that each vector I,II, “UL IV is 
of length m, It is also necessary to show that x = x, This is easily 
seen sincex@l=x and x @ l=x, therefore x =x (19 }) =x, 
Theorem 5.9 
Complementing one of the pair (I,/I) and one of the pair (III,IV) 


exchanges the Hamming vectors H(I,II) and H(IIl, IV). 


j 2385220202 n-] 
II=a ae 2oeeeeee ay 
Ill=b,b.,b,. oe ee b-] 
IV=b_b7¢ re ; bo 
eee | Nh e..... ae 
. HII, IV)=h)hoh3h, po ; he because of the parity check. 
Be H(1, I1)=H(1, I) =h)h,h,h, Soe hs 
a HI, WV )SH(II,1V) = hhh hea been, he hye 


by. Therefore H(I,II)=H(III,IV) =H(III, IV), 
6, and H(III,1IV) = H(, II) =H(1, I). 
ae ote aK 


ci 





Theorem 5,10 
Exchanging the (I, II) pair with the (IU,IV) pair exchanges their 
Hamming vectors. 


l, A change in the Hamming vector can only be caused by a bit 
change, 


2. Exchanging the (I,II) pair with the (III,IV) pair changes 
position but not bits, 


EKarlier in this chapter codes of length 40 generated from the quad 
were decomposed into kernels of length 10 by use of the Hamming 
vector. The conjecture was made at that point that perhaps codes of 
length 200 manufactured from kernels of length 10 could be decomposed 
into kernels of length 50. In order to prove this decomposition is not 
possible several more ee are necessary. These theorems will 
be proved before proceding to the decomposition problem, 

Theorem 5.11 

Under all possible general transformations of a complementary pair 
where n> 2 there are only four possible Hamming vectors, two of 
which are complements of the other two, 

IF (A,B) are a complementary pair of length n where n> 2, 


Zz. ee) = H interlaced with H where He are the odd bits 
of the Hamming vector and H_ are the even bits of the Hamming vector 


in reverse order, The symbol H=H. x ee will be used for A inter- 
laced with A in exactly the same manner that (I,II) and (UI,IV) are 


interlaced to form A,B. Also h. = a.@b, as before, 


mice oo 


2° Sane oe 


Ob =H Wow... My, 


oe, loet H,=(1@III)=a, %b), a.0b.. one Ob 


n=15 en | 


-.a58b 


H, =(i SIV)=a Sb , a. 39d _ x 


H,=(1@ IV)=a,6b_,a,0b, 5.....a, | 


H,=(1111I)=a ®b),a,_,8b,.....a,0b, | =K kA 


fave 





4. uv. =u.Ov. bye [Theorem tccuk. 


therefore 


(T@1TI)=H, (IIOTV )=H, (101 )=H, (11@ TIN)=H, . 
5. If h,=u,9v., then h,@1=h =u.dv.01=u.0v.=u.0Vv 


therefore 


rie = HU, V)=(U,V) 


and H, =(1@1II1) = (I@III) H, =(1@1V) = (I@TV) 
H,=(II@IV) = (II®IV) H,=(TI@ITI) = (II@TII). 
6. eee) = h hob, Weal 
Bee = hh | lee. Bae 
q But hea 2-7 n/2 roe by Theorem 5,6, 
8 H *H, =h) hy...----byyp hoyo..- hyhy 
H, *H, = Ei hoseeeee oh, Aare 22h 
9. Therefore H, H, = H, * Hy) = H, eH 1° 
@, ee El) = Ak, ‘eee Ki 
ae H, *H,=WM _)----- Ky. 
He. Ke /2- a ae 


3s = WT «< ty ae *x 
Me Therefore He, Hy A, k He, Ay H,. 


14, Table 5,1 is an exhaustive list of transformations ona 
complementary pair in the sequence quadruple form with the Hamming 
vector for each transformation, 


Field, Foo G 


se 
1 “Ha Hy H,*H 


1? Ha*H,» Hy*H 


3 and their complements are the 
only possibilities which exist, Therefore only four Hamming vectors 


are possible and by steps 9 and 13 they occur in pairs which are 


complements of one another. 


x 2 xe 


fs 





TRANSFORMATION 


HAMMING VECTOR 


COMMON VECTOR 


I (ot it IV) H(1)=H,*H, H(I) 
Cy (I II III IV) H(C =H," ‘A H(I) 
C, Cf lf til IV) HC, =Fi Jee (1) 
G (ie Til av) H(C) =H, *H> H(1I) 
ae ieeiiesitl 1) H(A, )=H,*H, H(I) 
A, (I It Ill Iv) H(A,)=H, *H, H(1) 
T) (lier iit. IV) H(T ,)=H,*H, H(T,) 
10 (I Il IV IID H(T,)= 3" ly H(T,) 
T (Ci ewes LT) H(T) =H, ir H(I) 
Zi (ier Le IV) H(Z) =H re ; H(T,) 
Y (II I Ill IV) H(Y) =H,* a H(T,) 
Xx (II I III IV) H(X) = H,*H, et) 
W (ieee tT) 1V) H(W) =H,* 1 H(T;) 
V (I II IV III) H(V) =H,*H, H(T,) 
U (tent sy f11) : H(U) ae Jel H(T,) 
S (I II IV III) H(S) =E Bel Bie 
R (I Ir IV Ip H(R) =H,*H, H(T)) 
Q (ity, oil) H(Q) =H,*H, H(1I) 
P (II I IV III) H(P) =H *H, H(I) 
O (ieeeti LV) H(O) =H, *H, Bhi} 
N (ote nil 1¥) H(N) = ease Et 
M (emery ~ TET) H(M)= H,*H, Ee 
Ib (lien TV) H(L) =H,*H, H(T,) 
K (I Il IV Iti) H(K) =H,*H, eee | 
J (II I IV_ It H(J) =H,*H, H(1) 
H (in iy 111) H(H) =H,*H, Eig) 
G (ie iy 111) H(G) =H,*H, H( 1) 
F Cieetho arr 1V) H(F) = H,*H, H(I) 
D (Cie Tiley) H(D) =H, *H, H(1I) 
T (II I IV III) H(1) =H,*H, H(1) 
B (Time LL) H(B)=H,*H, H(1) 
D oe Tih Bean) H(O) =H, *H, H(I) 


The E transformation just permutes the first two and last two pairs 

of (I II III IV) and since modulo 2 addition is commutative no change 
would take place inthe H vector. The last column identifies all 
Hamming vectors in terms of two vectors and their complements. 


TAB LE 
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Theorem 5.12 

If a complementary pair (A,B) of length n is operated upon bya 
complementary pair (C,D) of length m toform codes U,, Uz by the 
time sequence exponential method, or to formcodes V,, V2by the inter- 
lace exponential method and either U, or Ug is time reversed, symbo- 
lized by U orif Vi or V2 is time reversed symbolized by V then 
H(U) = H(V) = H(A’, B* yom where the (A? ,B') pair indicates some 
general transformation of the original (A,B) pair. 


A ® B = H(M). where B indicates the time reverse of 
the original B 
in A ® B = H(H). 

B ® A = A 2365 =H, ). 
4. B® A= A ® B=H(L). 


os Referring to Table 5,1, 


H(M)=H(T)), 


H(H)=H(T}), 
H(L)=H(T}). 
G c e d 
OF Ui a ae ck ca Lvaweuiee tine Bare 
— d = ll dy Cn cy 
Wee ae ek os A.» Bae. eee 
and 
C d E Cc d 
ee rs) m 
on Cn os | dy a) 
Ve=A B Pe eo EX Bas 


Time reversing one of the pair in each, for example, the use 


of Te yields 


ae 2 Ca oe 
Ui,=A / 2 eee oe ° A aR Pere 

C C C d d 
te= Bo BO B igs 

and Cy dy C5 Cn dn 

VizA A 6. eee A --B 

C d, c C d 
[a= 8) A Bee Big as 








le H(U) and H(V) consist of clusters of (a @ B 4) and 
d. d. 
clusters of (B° @ A ws 
Si The Hamming vectors for these clusters have the following 
four possible forms: 
A®B,A®B, B@A, andB @A, 
These four possible vectors were seen in steps 1,2,3,4,5 to all be 
equal to H(T,). If U; and V, had been time reversed, using the 
operation T ,, the result would have been H (Ta, 
EO). Since all the clusters are identical and there is a cluster for 
each of the m bits in C andeachof th m bits in D, 
H(U) = HV) = PHC rile for Uz or Ve time TWEE, BE 
imecene ral ; 
Vaan) 1 6 eld 
H(U) = H(V) = H(A’, B*) where (A ,B ) signifies some transfor- 
mation of (A ,B). 
*¢ x * 
Theorem 5.13 
The kernel of length 2 when operated on by a complementary pair 
by either the time sequence exponential or interlace exponential method 
can always be decomposed at least once by the interlace method, 
ie If one of the code pair formed by either exponential method 
is time reversed, H(U) =H(V) = H(A’, B?) ém by Theorem 5, 14, 
Zr There are only two possible Hamming vectors for the 
quad, Ol and 10. 


oe imeretore H{(U-) = Hivw)= fon or (10)*™ both of which 


are decomposable by the interlace method. 


sc 6 


* 
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Theorem 5,14 


A complementary pair (U,, Uz) formed by the time sequence 


exponential method cannot be decomposed by the time sequence method. 


Nh 


Let (U,, U2) be a complementary pair of length 2mn formed 
from the complementary pair (A,B) of length n operated 
upon exponentially by the complementary pair (C,D) of 
length m, formed in the same manner as in Theorem 5, 3, 

yn 


7 n n es mer | 
H(U)=(cy ®@ dm) (c,@d_ , -(c_ 6d, ) (d,@c oat} -ee(d_ Oc)) ’ 


The first mn bitsin H(U) are m clusters composed of n 
zeros Or n ones dependent upon the bits of H(C,D). 
Byiineorem 5,1 H=0'1™" or 17" 0" for decemipoeiaen 
of a code of length, 2mn by the time sequence method, 

This requires the distance of (C,D) be zero or m, but by 
Theorem 4.1 D(C,D) = m/2, therefore H(U) is not decom- 
posable by the time sequence method. 

ie onmor Joes time reversed sii(U)) = H(A), B1y*™ by 
ihieore ni Sec: 

This is incompatible with the form H=0°"71"™ since it would 
require the first m/2 clusters of H(A', B+) to have a distance 
of zero and the second m/Z2 clusters of H(A', B’) to havea 


distance of n, Both are impossible by Theorem 4,1 which 


requires a distance of n/2, 


ie * BAS 


Theorem 5,15 


A complementary code pair (V,1,V2) formed by the interlace ex- 


potential method (using the same definitions and symbols as Theorem 5, 4) 


can be decomposed by the time sequence method only if H(C,D)= 
gm/2, mi 2 a ym/e gm/2. 


fey 





2n 2n 2n 
Hey) = (cece! (c,0d__1) oe (c_9 d,) by Theorem 5, 4 


and equation 5, 5, 
The clusters of zeros and ones in this Hamming vector are in 


lengths of 2n and are dependent in value upon the bits of H(C,D). 


H=07 717" is the requirement for decomposition of this code 


by the time sequence method. 


m/2,m/2 


iia. bD)=0 or its complement, it will satisfy step 3 


since (o2nym/2 2nym/2_gmnymn_pyy) and it also satisfies 


Theorem 4,1. 


If Vi or Ve is time reversed H(V) =H(A?, B+)*™ by 


Theorem 5.12, 

This is incompatible with the form H=07""1™" since it would 
require the first m/2 clusters of H(A?,B’ ) to have a distance 
of zero, and the secondm/2 clusters to have a distance of n. 
Both are impossible by Theorem 4.1 since it requires a 
distance of n/2Z. 

Step 4 showed that H(C, D)zo™/4;M2_ This satisfies the re- 
quirement for the time sequence decomposition by inspection. 


Any other arrangement of H(C,D) would fail to allow decom- 


position, 


Theorem 5,16 


Any complementary sequence pair (Vi, Vz) formed by the interlace 


exponential method as in Theorem 5.4 (using the same definitions and 


symbols), the quad as the (C,D) pair, is time sequence decomposable, 


lips 
Ze 


The quad has only two possible Hamming vectors, Ol and 10, 


m/2)m/2 


This satisfies the form H=0 where m=2, the quad 


length. 
Vi and V2 are therefore time sequence decomposable by 
Theorem 5,15. 
oe x aK 
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Theorem 5,17 


A complementary pair formed by either exponential method is de- 


composable by the interlace method if and only if H(A, B4)=(01)"/4 on 


ao)/2, 


All definitions and symbols are assumed to be the same as 
those in W@iieorems 5, 3,5.4 and 5,12. 
i By Theorem 5,2 these codes are decomposable by the inter- 


a a aoe 


lace method only if they are of the form H=(01 

tag H(U) and H(V) are clusters of at least n ones or n zeros 
nm 
and therefore cannot have the form (01) = 


_— 


ae H(U) and H(V)= H(A?, B by Theorem 5,12, 


Hy 2 
Cc 


4, If and only if H(A‘*,B') has the form (o1)n/4 or (10) an 


my 2m _ (Oi orusing 10 instead 


H(U) or H(V)= [(0]) 
On a0 ie (Ones This satisfies the requirement of Theorem 
5,2 that the Hamming vector be of the form (01)" or (10)*, 
x 2 oe 
A brief description of the method used to prove that codes of length 
200 formed from kernels of length 10, cannot be decomposed into ker- 
nels of length 50, will now be given as a guide to the actual proof of the 
theorem. The method usedis to exhaust all possible methods of de- 
composition by examining all possible Hamming vectors for the (A,B) 


and (C,D) pairs, Table 5,2 lists all possible Hamming vectors for 


Kernels of length 10. 


Kernel 1 Kernel 2 

A=1000000110 A=1001101111 

B=1001010001 B=1100001010 

The 4 possible Hamming vectors, The 4 possible Hamming vectors. 

0001010111 0101100101 
1110101000 1010011010 
0000101111 1100101100 
1111010000 0011010011 


All possible Hamming vectors of kernels of length 10, 
TABLE >. Zz 





The two methods of generating codes of length 200 from kernels 
of length 10 are the time sequence exponential method and the interlace 


exponential method given respectively by 


Cc C Cc d 
en a “eae oR DS B10 
d d d. ¢ C (er0)) 
ee A A Be 
Cc d Cc d Cc d 
Wot MEI Bs, ee A Bee oo 
d C d C d Cc : 
Viele tO) eee i bee? 


where (A,B) and (C,D) are kernels of length 10. 

The Hamming vectors of these codes of length 200 are anti-sym- 
metric about their centers andare inclusters of 10 zeros or 10 ones 
for H(U), andin clusters of 20 zeros or 20 ones for the H(V) vector. 
Time reversing either U,; or Uz and similarly either V, or Vos , 
gives 20 clusters of Hamming vectors of the forms listed in Table 5,2, 
since each of the 20 A codes is now matched with some transform of B. 

The known possible ways for kernels of length 50 to generate com- 
posite codes of length 20 are: 

IF. Either interlace or time sequence kernels of length 50 to 

form composite codes of length 100 and then either inter- 


lace or time sequence these n=100 codes to form composite 
codes of length 200. 


rig Take kernels of length 50 and either time sequence expo- 
nential or interlace exponential with the quad to form codes 
of length 200, 


oe Take kernels of length 2 and either time sequence exponen- 
tial or interlace exponential with kernels of length 50, 


An examination of all possible Hamming vectors of length 200 from 
each of these methods will show the incompatability of these Hamming 
vectors with those generated by the kernels of length 10, 

Theorem 5,18 

No codes of length 200 formed from kernels of length 10 canbe 


decomposed into kernels of length 50 by standard methods of decomposi- 


tion. 80 





Given U, Uz and Vi Vz both complementary pairs of length 
formed by the time sequence exponential and interlace expo- 
nential methods as given in equations 5.6 and 5.7 respectively. 
The time sequence method of decomposition is impossible for 
each of the following Hamming vectors for the following reasons. 
a. H(U) can never be decomposed by time sequence. 


Theorem 5,14, 


lo H(V) requires a form Nei by Theorem 5,15, 


Table 5,2 contains no vector of this form. 
Gx H(V) and H({U) can never be decomposed by time 


sequence, Theorems 5, 14,15, 
The interlace sequence method of decomposition is impossible 
for each of the vector forms for the following reasons. 


a. H(U) and H(V) can never be decomposed by interlace, 


Theorem 5,17, 


b H(V) and H(U) are decomposable by interlace only if 
es, ee z 


H(A?, B’) is of the form (01)> or (10) by Theorem 5,17, 
Table 5.2 contains no such vectors, 
Assume that one of the Hamming vectors H(U) or H(V) hasa 
second exponential form, first with the kernel of length 2 
operated on exponentially by the kernel of length 50, and 
secondly by the kernel of length 50 operated on exponentially 
by the quad, Let H’* be the designator for the Hamming 


vectors of the 2 and 50 combinations, 


as Consider first the quad acted on exponentially by the 
kernel of length 50, HY U)= H(V) = H*(A',B’*). But since 
A’,B is the quad, by Theorem 5,13, it is decomposable 
by the interlace scheme. This is not possible since 
step 3 above exhausted all possible interlace decompo- 


sitions, 


bie Next consider the kernel of length 50 acted on exponen- 


tially by the quad. Considering first the exponential 
8 1 





interlace scheme, by Theorem 5,15 if Hi (GD )=0 loom aon 
which it does for the quad, Therefore the code is time se- 
quence decomposable, However, in step 2 all of these time 
sequence possibilities were exhausted ard therefore this de- 
composition is not possible, H*(U) would have a form con- 
taining either 50 zeros or 50 ones inclusters, For the 
original 10 length codes to form these would require a 
0>1° Hamming vector which is not listed in Table 5,2, 

This exhausts all possible Hamming vectors that might lead to de- 
composition, therefore codes of length 200 formed from kernels of 
length 10 cannot be decomposed into kernels of length 50 by standard 
methods, 


ye oe aK 
Table 5.3 shows the possible Hamming vectors for codes of 


length 26. 
A=01001101111010111100111010 
B=10110010000111111100111010 
The 4 possible Hamming vectors, 
11111111111101000000000000 
OOODODODNDOOOOOIOLIIIILIIIANI 
11101110111010100010001000 
00010001000101011101110111 


All possible Hamming vectors of the kernel 
of length 26, 


table Seo 
An examination of Table 5,3 which is exhaustive for the Hamming 
vectors of length 26 and Table 5,2 for 10, shows that the next two 
theorems can be proved by identical methods to those used in Theorem 
5.18. Since these proofs exactly parallel that of Theorem 5,18 they 
will not be given, 


Theorem 5.19 


Complementary pairs of length 520 formed from the time sequence 
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exponential or interlace exponential from kernels of length 26 and 10 
cannot be decomposed by standard methods into kernels of length 130. 
xe aK oS 
Theorem 5, 20 
Complementary pairs of length 1352 formed from the time sequence 
exponential or interlace exponential method from kernels of length 26 
cannot be decomposed by standard methods into kernels of length 338, 


se oe oc 


The study of Hamming vectors thus far has disclosed that a pattern 


Pao coon! © nt 


or their complements in either the (A,B) 
pair or (C,D) pair, is required in order that a code formed from one 
kernel length by an exponential method be decomposed into a pair from 
a different kernel length, The quad is only kernel which thus far satis- 
fies this criteria, However, the other kernels come Mee as a check 


shows that pele and (01)"10(01)° are among the Hamming vectors of 


anon S is a Hamming vector of length 26. 


length 10, and 0 

If the quad is involved in the formation of composite codes it makes 
their decomposition more likely. This is shownin Theorem 5,12 and 
Theorem 5,16, Appendix II demonstrates this fact quite clearly since 
it contains all codes of length 16, less some operational redundancies, 
Each of the code pairs and its time inverse is decomposable by some 
method into shorter length codes. Appendix III] lists all the pairs of 
length 20, less some operational redundancies, Each code pair or its 
time inverse can be broken into its originating pair through the use of 
Hamming vectors, but no code and its time inverse can both be decom- 
posed by the standard Hamming vector methods. 


The writer would like to end this chapter with a conjecture that 


kernels of length 50 do not exist, This conjecture is based upon three 
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bits of evidence, the proof of Theorem 5.18 being the first, The second 
piece of evidence is based upon the fact that 50= Tae Sacro From this 
one might assume that codes of 200 formed from 10's might cover just 
one set of the ones determined by equation 2,5, This proved not to be 
true however, since the number of ones in both the exponential formations 
of codes of length 200 from kernels of length 10, if they were decompo- 
sable into kernels of length 50, would exhaust all possible unordered 
pairs of ones for n=50 as determined by equation 2.5, The last bit of 
evidence to base this conjecture upon is that the code length is not twice 

a prime number. The only previous kernel which might have existed, but 


was not twice a prime number, was n=18, and it did not exist, Similar 


reasoning extends this conjecture to kernels of length 130 and 338, 


84 





CHAPTER VI 
SUPPLEMENTARY AND CYCLIC COMPLEMENTARY CODES 

This chapter is concerned with two classes of codes, The first of 
these classes consists of quadruples of sequences with the property that 
the total number of likes at each spacing equals the total number of un- 
likes at the same spacing. These are called supplementary codes. C The 
second class of codes, which will be defined later, consists of the cyclic 
complementary codes, Both of these code types have complementary 
sequence pairs as a possible subset; this allows the use of their prop- 
erties as screens inthe search for new kernels. 

Supplementary codes will be discussed first, since the nee search 
involved the use of this property as a screen, while the cyclic comple- 
mentary property was not used until the search for kernels of length 34, 
Only one theorem will be proved for each of these types of codes, since 
these two theorems were the only ones applied in the actual search for 
new kernels. 


ibe UpUpUgeseeeee UL 


Non De ee 2 eee 


n 
Ww WwW Ww eo @ 8 8 @ WwW 

aed, => n 
Kp XoXgeeeece- x are sequences of zeros and ones and 


satisfy the constraint of equation 6.1, they form a quadruple which is 


supplementary. 
i= j 
y cee asi 7°? + DD aes ot) + Wee + eer = 
i=l 
j 
AO geatig)! FY, OMe ay Be eae ee (6. 1) 


=| 


for all j, 1<j<n-1. Note that there are 2j like pairs and 2j unlike 


pairs at a spacing of n-j. 
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For simplicity in the proof, we use the symbolism that was used in 


Chapter 2, 
J 
let Li; - u. 7] nest @ | tor all’ j,i <7 sen 1 
j 
“aes ti -j 


A restatement of equation 6,1 inthis symbolism is 


Ly,t+Lyt+hly+lLy =Uy+tUy + Uy + Uy (6, 2) 


Theorem 6.1 

Any complementary sequence pair, (A,B), written in standard 
(I II III IV) form has the property that I,II, III and IV area 
supplementary quadruple, 

it LA= B and U 


oe Considering just the even values of j rather than all values, 


iG by the definition of complementary. 


A “B 


L,=L, 7 


UT ats Up 


Derby t yy» 


IJ’ 


U + U, 


BOW ie 


a L, 1 Li =Unsy 15 Uy: 


4. LiythyyF¥ + Ur 


2 Adding the equations in steps 3 and 4 gives 
L, + Ly + Ly, + Liy = U; + Un 7 Un + Un : 
6. This is the same as equation 6.2 and I, II, III, IV formed 


in the sequence quadruple form from a complementary pair 
are supplementary. 
Since the supplementary property applies only to the even spacings 
in a code pair, the converse, that all interlaced supplementary quad- 


ruples are complementary, is not true, 


sk oc xk 
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The cyclic complementary property is useful in itself for a communi- 
cations or telemetering system; however, at this time our concern is 
for its property as a necessary condition for a code pair to be complemen- 
tary. Acyclic sequence or code, as the name suggests, 1s a never 
ending sequence of zeros and ones which has a period of n bits, A cyclic 
complementary sequence pair, is a pair of cyclic codes, each of period 
n, where the number of likes of one sequence equals the number of un- 
likes of the other sequence for n possible matches for all spacings 


from 1 to n-l, or stated interms of the bits of the A and B codes 


n sles al 
oy = a a. ae a er eee he ] i=) < no (6. 3) 


where a ..=a. and b 
nt1 1 n 


Theorem 6.2 


4 =b. since the periodis n, 


A complementary sequence pair (A,B) of length n, if writtenina 
cyclic fashion, is always a cyclic complementary pair of period n, 


das Since (A,B) are a complementary pair they satisfy 
equation 2,3. 
. i: a 
i ve a. Smee ae b.® Lect 1 for j, l=j=n-1. 
ae Expanding for fy gives 


a,®a <=b,8b 61 whichis equalto a $a, =b 6b,6 1, 
n n n n l 


] ] ] 


oy Expanding for fi] gives 
(a,® as) + (a,® a2) + a.® ay) eae (a i® a_,)=(b, @b,@1) + 
(b,@ b,@ 1) + (b,0 b, 1) +......+(b, 9 6,8 1). 

4. Adding the equations of step 2 to step 3 gives Ci. 

2) Similarly expanding f, gives 
(a6 a_i) + (a5@ a) = (b,@ b,_19 1) + (b5@ b® 1) 

which is equal to 
(a__ 8 ay) + (a 6 a5) = (b__,9 b,® 1) + (b_® bo OF le 
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On Expanding a gives 


(a, %a,) r (a,0a,) + (a,0a liste (a, _2%a_) = (b 6b,6 ye 


5 


(b,@b,@1) + (b,@b,@1) +....+ (b 


‘a Adding the equations of step 5 to step 6 gives C,. 


2b, 61) 


8. Continuing this same procedure through all possible values of 


jJ would show that oe = - of Therefore the cyclic com- 


n-Jj° 
plementary constraint equations are based on the sum of two 
restricted portions of the complementary constraint equations, 
Therefore all complementary pairs have to be cyclically 
complementary. The converse is not necessary true since 

o is a sum, and ina sum the addend and augend are not 
unique, 


uK KK * 


Since ws = = a fete each constraint of the complementary property 
is used twice as j varies from 1 to n-1; therefore Oo. is symmetric 
and is centered at Ci /2: In the application of the cyclic complementary 
property in Chapter 8, half of the characteristic cyclic number is de- 
leted due to this symmetry property. 

Both of the screens described in this chapter were suggested by 


ce ee 


Dr, Golay; however, to the best of my knowledge the theorems 


and their proofs have not appeared in print, 
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AN EXHAUSTIVE SEARCH FOR KERNELS OF LENGTH 26 
CHAPTER VII 


In his paper Golay emphasized the need for an exhaustive search for 
complementary code pairs of length 26. This chapter describes sucha 
search, The purpose of this search was to determine if there were any 
kernels of length 26 and if any of these exist, the number of such 
kernels, It would prove very difficult to surmise the number of kernels 
in a code of this length predicated on the known kernels of other lengths, 
samniee n=Z2 had one kernel, n=10 had two, and n=!18 had zeron Anex_ 
haustive search for kernels of length 26 might throw some light ona 
possible general method for finding kernels of longer lengths, or might 
disclose some sort of pattern showing the distribution of ke ele among 
the possible code lengths. 

The only feasible method to accomplish this search was with a high 
speed digital computer. This chapter is concerned with the computer 
program andthe results of the exhaustive search for kernels of length 26, 

At first glance it would appear that there are 224 possibilities to 
be screened for an exhaustive search for all possible kernels of length 
26. Although this is true, the application of some of the theorems de- 
veloped in the earlier chapters immediately eliminates from considera- 
tion large blocks of the joe possibilities which are either redundant or 
impossible. This is a very necessary procedure since even with the 
highest speed computers of today it would be impossible to investigate 
this number of possibilities in a life time. 

The program was developed for the CDC 1604; this computer ac- 
complishes approximately 200,000 operations per second, This speed 
coupled with the powerful screens used in programming made possible a 
reduction of the computer run-time to approximately 75 hours, A brief 


outline of the various screens used and their reduction factor will now be 
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given before continuing in more detail with the actual programming 
techniques used, The word code or code pair as used in this chapter 
and also in the next chapter will be understood as possible complemen- 
tary code pair, 

The first screen utilized was the number of ones which must appear 
in each of the code pairs. This was given by n= (aera + een 
as derived in chapter 2, where n is the length of the code, p is the 
weight of the A code and q is the weight of the B code. The possible 
solutions for n=26 were the unordered weight pairs (16,15), (16,11), 
(11,10), (15,10). The pair (p=16, q=15) was arbitrarily chosen to be 


used in the program. This reduced the total number of possible code 


Be is: fas oy 13 


od 


patrs trom 2 35x10 down Pee Tir ~ 4x10 ~ . This 
was a reduction bya factor of 200, which brought the life time search 
down toa little less thana year, but this amount of time was still not 
feasible for a computer search. Golay remarked ina footnote thata 
complementary pair might be thought of as being composed of two inter- 
laced half length codes, in this case n=13, 7 Since it was obvious the 
total number of ones in both of the interlaced pair must equal the number 
of ones in one of the complementary code pair, equation 2.5 was 
applied again, this time for n=13, The unordered weight possibilities 
for n=13 were (9, 7), (9, 6), (6,4), (7,4). In order to be compatible with 
the (16,15) pair for n=26, the (9,7) pair was selected for the code con- 
taining 16 ones and the (9, 6) pair was chosen for the other code which 
contained 15 ones. 


The basic problem had now been reduced to all possible combinations 


of the four thirteen bit sequence families, This was equal to 


ie 3. i. 13. . 
jee ea * agi x pr = 715x715x1716x1716 or approximately 
12 


15x10 possibilities. Although this reduced the problem by one more 
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order of magnitude it was still much too large for a practical search, 

Dr. Golay ina private communication with the author pointed out the 
possibility of using the supplementary characteristic as an additional 
screen, oe Theorem 6.1 gives the necessary and sufficient condition 
that complementary codes sequence quadruple form be supplementary 
sequence quadruples. The codes of length 13 were therefore categor- 
ized according to their first 3 and last 3 bits. The number of like 
pairs at spacing twelve, at spacing eleven, and at spacing ten was com- 
puted for each of the 13 bit sequences, For example, 0111100111011 
has 0 likes for spacing twelve, 1 like for spacing eleven, and 3 likes 
for spacing ten. The block number given to 0111100111011 was there- 
fore 013, Out of sixty four possible combinations only eighteen different 
block numbers were generated for each of the twS number sequences of 
all possibilities of thirteen bits with nine ones and with seven ones, An 
examination of Figure 7, 3 reveals that there are a widely varying number 
of members within each of the blocks, The supplementary characteristic 
states that the total number of like pairs for each spacing must equal the 
total number of unlike pairs for the same spacing, With four codes ata 
spacing of twelve there are 4 possibilities, therefore two of them must 
be likes; at a spacing of eleven there are eight possibilities, therefore 
four of them must be likes; and likewise at a spacing of ten there are 
twelve possibilities of which six must be likes to satisfy the supplemen- 
tary property. Block numbers were added for each possible combination 
to see if the total were 246, and in cases where the total was 246 those 
particular codes within these blocks were then sent to the next necessary 
but not sufficient condition for the codes to be complementary. If the 
total was not 246, the blocks of codes were rejected. The power of this 


screen is shown in the example where say two of the blocks had code 


e) I 





numbers 000, 000. The only possible combination to add to 246 
would be 123,123 for the other two codes. This cuts the total number 
of possibilities in this particular block combination down from 


6 
“ to only 64 possibilities. Although this admittedly 


1716x1716 = 2, 9x 
was the most extreme example, this screen reduced the number of 
possible codes checked to about 107 or further reduction of three 
orders of magnitude. 

The next necessary condition to be checked was the parity test, 
equation 2,4, All code quadruples which satisfied the parity test were 
then sent to the necessary and sufficient like pair, unlike pair check for 
each spacing. The like, unlike check subroutine was quite long and also 
involved word unpacking. This subroutine which was ignored in the 
calculation of run time increased the computer search time from a calcu- 
lated 70 hours to an actual runtime of the order of 75 hours. 

The half length codes in the program were interlaced in the same 
manner as was done in the operations group formulation, This method 
of combination made the parity check easy to calculate since the bit 
positions in the computer words were the same as the sequence posi- 
tions of the codes, Modulo 2 addition was directly applied to the code 
quadruple to forma vector of all ones, providing the parity test held. 

If the code vector was not all ones the code quadruple was rejected. 

It is to be noted that by selecting the ordered (16,15) pair for the 
number of ones in the A and B code, the operations Ci, Cz , and C 
were eliminated from possible redundant consideration, These were 
eliminated by their requirement for (10,15), (16,11), and (10,11) 
ordered pairs of ones respectively, based on equation 2,5, The inter- 
lacing of the half length codes with ordered pairs of (9,7) ones and 


(9,6) ones respectively eliminates from consideration the operations 
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Ti, Ts, T and A, since these would have required (7,9,) (9,6), 

(9, 7)(6,9), (7,9)(6,9) and (4,7)(4,6) ones respectively. The only basic 
transformation not deleted by controlling the weight 1n each of the code 
quadruples was Ag. This is seen to be true because if the pair 
(9,7)(9,6) has the operation Ag performed on it, the result is (9,6)(9, 7). 
If this code is now transformed by the trivial exchange operation, E, a 
code pair of the original form is generated, The scheme which was used 
to eliminate the Ag redundancy will be discussed late in this chapter. 

Figure 7.1 isa rough flow diagram for the program used in the 
exhaustive search for kernels of length 26. With modifications in data 
as necessary, this program was also used for exhaustive searches of 
composite codes of 16 and 20, 

A brief description of the blocks in Figure 7.1 will give the reader 
some familiarity with the programming philosophy used in this search, 
The actual programs as written for the CDC 1604 are given in Appendix V. 

The number generation subroutine is able to generate any number up 
to 48 bits in length (word size of the 1604) with any number of ones up 
to a maximum of 48, This is accomplished in the 1604 by inserting a 
number of the form 00000...Q00001111...111 into the accumulator, 
where the total length of the number is 48 bits and it contains p ones 
and has a code length of n, This number is then shifted left 48-n bits. 
This shifted result is then checked to see if it is negative, which indicates 
a one in the most significant bit position, A procedure is now set up to 
count a one or zero and then shift left one bit and repeat the process for 
the entire n bits. If p ones are counted the number is stored, if not it 
is rejected, A number one larger than its predecessor is inserted into 
the system and the process is repeated until the ones are all in the 
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Blocks 2,3,4 are used to compute a supplementary characteristic 
number for each of the 715 thirteen bit numbers with 9 ones and also 
for each of the 1716 thirteen bit code numbers with 7 ones; then in 
accordance with their supplementary numbers these are stored into 
group blocks, The supplementary number was derived by using a mask 
to expose for consideration only the first 3 and last 3 bits of each of 
the 13 bit code numbers, These masked numbers were then compared 
with an exhaustive list of all possible combinations of zeros and ones in 
these 6 bits. These 64 possible combinations were in an M to one 
correspondence with a list of 18 possible supplementary characteristic 
numbers (see Figure 7.2). This correspondence allowed the temporary 
attachment of this supplementary number to each of the codes, The code 
numbers were rearranged into blocks according to supplementary number 
and at the same time had their temporarily attached supplementary 
number deleted. This part of the program is purely for putting the data 
into a useable state and is used only once per run, The rest of the pro- 
gram is highly iterative. 

Consideration of the (I II IH] IV) form of the possible 26 bit 
codes, shows that each of the 18 blocks of I, II, Ill, and IV must be 
compared against each other indicating 18x1!8x18x18 ~105,000 possible 
block comparisons, This was cut down somewhat by use of a slightly 
different technique, This was accomplished by loading the computer 
accumulator with 246 and from this subtracting the first block number 
of I, then subtracting the first block number of III, and then subtrac- 
ting the first block number of II, Rather than subtracting the Block 
number of IV, the difference obtained was instead checked against all 
all possible block numbers of IV to see if it was listed. Since if the 


block number did exist only one could exist. This procedure cut the 
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possible number of block comparisons by a factor of eighteen and left 
only 5800 combinations to examine rather than 105,000. If the block 
did exist all the codes within this {I II III IV) block grouping were 
then giventhe parity test, If the IV block number did not exist II was 
stepped ahead one block number and the process repeated. After II 
had cycled throughits 18 block numbers, II1I was then stepped ahead 
one block number and the entire process repeated. Similarly when III 
had cycled its 18 block numbers, I was stepped ahead and the process 
repeated. To avoid the Ag operation redundancy, whenever I was 
stepped up one block number IIL was started from this same block 
number rather than at the first block number, This avoided ate bioek 
number combinations (I III),and (III I) both appearing except when 
dlock I was equalto block III. This block number cycling allowed a 
vantage breakpoint for partitioning the program into suitable size in- 
crements for computer run times, 

The parity check was made by adding(Modulo 2) a code from I 
to acode from III toa code from II, This resultant was then checked 
to see if it existed as a code inthe list fromthe IV block, If it did 
exist it was complemented and then sent to be unpacked for the necessary 
and sufficient like, unlike test, If it did not exist another set of codes 
was sent in and the process repeated until all possible code combinations 
in (I II ID) were exhausted. 

The methods used in the unpacking and like-unlike subroutines are 
quite straight forward and will not be amplified here, although it is to be 
noted that the like-unlike subroutine is not limited in the lengths of the 
codes it can test, up to the machine storage size, whereas all the other 
subroutines, since they use packed words, do have quite restrictive 
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iris program was first used to make an exhaustive search for 
codes of length 10 for test purposes since all kernels of this length 
were known, After this satisfactory checkout the exhaustive search for 
codes of length 26 was undertaken, 

During the first part of the computer runfor n=26 Dr. Golay in- 
formed the writer by personal correspondence that he had discovered a 
kernel of length 26, using a "by-hand" Pechnictieta at The exhaustive 
search revealed only a transformation of the code pair that Dr. Golay 
had discovered. Therefore it was proved that only one complementary 
kernel of length 26 existed and that it was (ignoring allowable trans- 
formations) 


A=01001101111010111100111010 
B=1O1L1LOOLOOOOLIILIILIOOLLIOIO , 
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CHAP TIC Retay vel 
A PARTIAL SEARCH FOR KERNELS OF LENGTH 34 


The next possible kernels for investigation after length 26 were 
those of length 34. As will be shown in this chapter an exhaustive 
search for kernels of length 34 was not feasible. However, an impor- 
tant subset of possible kernels was exhaustively searched and a scheme 
which can be adapted to the general search was programmed, No ker- 
nels were found in this partial search. To use the same method of 
attack as that used on n=26 did not seem too feasible, both because of 
the length of time involved in the search and because there would be no 
guide to searching, beyond random _ Picking of possibilities until all 
possibilities were exhausted, The length of time involved for an ex- 
haustive search of n=34 was estimated fairly reasonably by checking 
the number of possible codes of length 26 against the number of pos- 
sible codes of length 34, since the time for the 26 search was known 
with reasonable accuracy. For n=26 there were 715x715x)716xl1716~ 
A ssell(p code possibilities at one point inthe screening process. At 
the same point in the screening process for n=34 codes there were 
12376x123 76x19448x19448 = 5, gx19l° possibilities, which is a ratio of 
4x1 o# to one as compared to n=26. A time estimate from the 26 
length code to the 34 length code was therefore fete 10 =3510° hours, 
This length of time was of course not practical for a search of an ex- 
haustive nature, 

The use of the cyclic complementary property as a screen was 
pointed out by Dr. Golayas highly practical because some reasonable 


Ie 


guess could be made as to which cyclic sets might contain a kernel. 


This would lead to a hybrid ''by-hand'' and computer search, It would 
also be easier to document the areas searched by the cyclic method as 


compared to the supplementary method and would therefore avoid the 
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duplication of effort by future n=34 length searchers. As will be dis- 
cussed at some length, there is in addition an even more important 
reason leading to the decision to substitute the cyclic complementary 
property for the supplementary property as a screen. 

The generation of cyclic sets which are to be utilized in the search 
for new kernels, is accomplished in a somewhat unusual fashion. This 
method of classifying, puts within the same set all possible numbers 
which are cyclic permutations of one another, as is to be expected; 
however, there are also contained within the same set all numbers 
which are formed by the removal of every other bit in a cyclic fashion 
until all bits are used, As an example the five bit number 11010 will 
be used to show all the numbers contained within its set, First, 


selecting every other bit, starting with the first bit, and repeating on 


the result gives: le 11010 
ap 10011 
De 10101 
4. 11100 
Be 11010 


and second, by permuting the bits of each of these in the normal cyclic 


permutation manner yields 


1. 11010, OO lL. e0 FO Ooo OA 
ip LOGI, LlOCi arti, O11 0. SCOT t: 


It is not necessary to generate the cyclic permutations from 3 and 4 
Since they are contained within 1 and 2, 

Counting cyclically the number of likes for each spacing gives the 
code number 1331 for 1] andthe code number 3113 for 2, These 
code numbers start with a spacing of 1 and include up through a spacing 
of n-1l or four in this particular example, The code number 3113 
means that at a spacing of one there were 3 likes in 10011 counted 


cyclically, at a spacing of two there was 1 like, at a spacing of three 
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there was 1 like, and ata spacing of four there were 3 likes. The 
example just givenis somewhat trivial since all binary numbers of 
length 5 with 3 ones fall into the same set, However, this is not the 
case for n=26 as there are 7 sets for codes of length 13 with 9 
ones and 76 sets for codes of length 13 with 7 ones, 

An examination of the kernel of length 26, shows that its sequence 
quadruple of codes came not from four different cyclic sets as one 
would expect but came instead from just two cyclic sets, 

Using the kernel form 


A=01001101111010111100111010 
B=10110010000111111100111010 one can decompose this into 


LOOP CTIIT O11 
II=0010110011101 
Piet 10 1 OCu ro 1AM 
IV= 0010111100010 


Using I to form its cyclic set gives: 


~ 9010111110111 
emo LLL LOOT 01 
ree OSI OE allan LEO LG, 
OO 1 OOMMET I 1 1 
meeOOlOl) blekot)! 


— Pb Ww PP 


Now take 3 and write it twice, this gives 01110111(1101001110111) 
11010. Picking off from the ninth to the twentieth second bit gives a 
code which is the sameas III above. Therefore I and III are 
members of the same cyclic set. Similarly take II and form its 
eyelic set; this pives: 


0010110011101 
~ OLLOLILOOIL0OL0 
0111000101011 
s OFOOCOOITITO)T1 
~  OOOTLTITOe1 01 


om BP WW fo KH 
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6 0011011011100 
7 0101110011010 
8. 0010100111011 
9. 0110101000111 
ieOrS Ola O11 100001 
ies Orel OOl)1 100 
rae OCOTLIOIITOLIO 
1, 0010110011101 


To make a comparison between II and IV the complement of IV will 
have to be used in order that it will have 7? ones as does II, 

IV= 1101000011101 

Writing 4 fromthis list of 12 twice gives 01000011101(1101000011101)11. 
Picking off from the twelfth to the twenty fifth bit gives IV, which in- 
dicates that II and IV are from the same cyclic set, 

A check of kernels of length 10 to see if their quadruples came 
from cyclic sets in pairs, as the 26 kernel did, proved that both ker- 
nels of length 10 didin fact come in pairs. This was meaningless 
however, since as was seen in the example of length five with three 
ones, there was only one cyclic set possible for each pair. 

In the hope that if the kernel existed inthe n=34 case, it would 
also be formed with I, III] from one cyclic set and II, IV from another 
cyclic set, the program ae modified, This change along with giving an 
exhaustive search for code possibilities taken two at a time from cyclic 
sets would also give an estimate of the computer run time for the ex- 
haustive search in general using the cyclic complementary property. 
This was the primary reason for the change in the program. 

Eight octal digits were used to designate the cyclic complementary 
number for each of the 12,376 possible half codes of length 17 with 
11 ones and for each of the 19,448 possible half codes of length 17 
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with 10 ones. Only eight digits were necessary as each cyclic count 
of spacing i actually counts both the i spacing and the ntl-i spacing 
simultaneously, This was demonstrated by step 8 in the proof of 
Theorem 6.2. The example of the code of length five given earlier in 
this chapter showed this symmetric property, as the code numbers were 

1331 and 3113, The first two digits of these cyclic complementary 
designators contain all the information available. 

The revision made to the computer program was not too extensive 
as only blocks 2,3,4,5 in Figure 7.1 were modified to handle the 
@yelic complementary property rather than the supplementary property. 
Of course all blocks were modified to handle codes of length 34 rather 
than codes of length 26. . 

One difficulty encountered was the lack of computer memory for the 
large blocks of data, and it was necessary to store the cyclic sets on 
magnetic tape rather than in the main memory core as the 26 case was 
handled, The cyclic sets were then called in two ata time for checking. 
This slowed the search down a little, but by no more than 1% of the 
total time. Figure 8.1 is a block diagram of the process of obtaining 
the cyclic complementary numbers and it will be worthwhile to remember 
Hoe memory size restriction when reviewing the procedure used, 

A general description of Figure 8.1 will show the differences in 
the method used in the search for n=34 kernels as compared to that 
used for n=26 kernels, The procedure inthe n=34 data formation was 
identical with that of n=26 as far as the initial possible number genera- 
tion was concerned. (Block 1 in Figures 7,1 and 8.1). There was 
an immediate departure from the old method in blocks 2,3, and 4 
where one of the 17 bit numbers with 11 ones was shifted left 17 bits 


and then added to itself, This formeda 34 bit number which was 
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actually the 17 bit number written twice in time sequence. A mask 
was then put over the last 17 bits and this masked number was used, 
with an equality searchof the 12,376 generated 17 bit numbers, to 
identify this particular number inthe list, After the number equal 
to the reference number was located in the generated list it was tagged 
a one, indicating cyclic set one. The whole 34 bit number was then 
shifted right one bit and a mask used to again pick off the last 17 bits. 
The search through the list and tagging with the group number was re- 
peated. After cycling through all 17 possible cyclic codes and tagging 
each, the original 34 bit number was now repeatedly long right shifted 
one bit and then right shifted one bit. A long right shift saves the bit 
which is pushed off the A register and stores itinthe Q register, 
while the right shift just pushes the extra bit off the end of the A reg- 
ister, Therefore every other bit of the 34 code was saved in Q and 
this formed a new 17 bit code within the same complementary cyclic 
set. This code was then formed into a double length code and the list 
of all possible 17 bit numbers with 11 ones was searched and tagged 
with a one and the process which was first described to generate cyclic 
permutations was repeated. This process of first forming 17 cyclic 
permutations of a code and then taking every other bit, was continued 
until the possible code which was originally operated upon reappeared, 
Then cyclic set 2 was started and tagged accordingly, and so on 
through 98 sets. The 19,448 codes of length 17 with 7 ones were 
classified in exactly the same manner and formed 150 cyclic comple- 
mentary possibility sets, 

After each series of sets were formed a mask was used to search 
the set tag numbers and pick off each set for transfer to magnetic tape. 


To insure that the pick-off was exhaustive, a zero was inserted on the 


list in place of each removed tagged number. 
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The sets were now ready to be called back into the computer for use 
in the search. Normally only one set was called in at a time because 
the pair fromthe I, IH set was checked against allthe II, IV sets 
before a new I, III set was called into the computer, The I, III sets 
contained the 11 weight codes while the II, IV set contained the 7 
weight codes. Whena set was called into the computer each cyclic 
permutation subset of 17 possible codes was givena cyclic comple- 
mentary number to categorize the entire subset, Each group contained 
either 4 or 8 of these cyclic subsets, These subset numbers were 
formed by counting the likes forthe I, III codes and the unlikes for the 
II, IV codes for the eight unique spacings. The codes within the subsets 
were checked against each other only if the likes equaled the unlikes of 
the cyclic complementary subset numbers for all spacing, or Ly + 
ies 0 ee 

Figure 8,1 is joined to Figure 7.1 just beyond block 5 and takes 
the place of blocks 1,2,3,4,5 in Figure 7,1. The same feedback paths 
shown in Figure 7,1 to insure the checking of all codes within a set 
and to feed in a new set after a check has been completed, are still in 
operation for the same purpose, Appendix VI contains the computer 
program for Figure 8.1. 

The total run time for the n=34 search, taking I, III from the 
same cyclic set and II, IV from another cyclic set was approximately 
15 hours. No codes were found. An extrapolation from this run time 
to the exhaustive run time for n=34 gives 15x150x98 ~2, 2410" hours. 
which although one order of magnitude less than the supplementary 
method, is still not reasonable, 

A second partial search for n=34 kernels was attempted after 


observing that the kernel of length 10 
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A=0101000011 
B=00001001 10 


with its center quad removed yields 


A*=01010011 
B+=000001 10 


which is a complementary pair of length 8, Breaking A‘*,B? into its 
(I II III IV) configuration gives 


I=0001 
II=1011 
IiI=0001 
IV=0100 


where (I, II) and (III, IV) are both complementary pairs, 

An attempt was made to draw an exact parallel from this observa- 
tion to codes of length 34, ‘One possible solution to equation 2.5 which 
determines the number of ones necessary in each of the kernels of 
length 34 is (21,18). If this same equation is applied to codes of 
length 17, two of the possible solutions are (11,10) and (11,7). The 
(11,10) pair when interlaced will total 2! ones, and similarly the 
(11,7) pair when interlaced will total 18 ones, Whereas the kernel of 
length 10 with the quad removed had a sequence quadruple form with 
codes of length 4, the 34 code with the quad removed would have a 
sequence quadruple form with codes of length 16. Appendix II gives 
a list of all possible codes of length 16 with weight 10 in the A code 
and weight 6 inthe B code. This list does not include the results of 
the time reversal transformations which have to be included since weight 
is invariant under T, , T2and T. It was observed that if these codes 
of length 16 were interlaced to form codes of length 32 and then sep- 
erated in the middle and the quad A=l1, B=10 inserted, all conditions 


as far as the number of ones required in both interlace codes of 17 to 
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form codes of length 34 would be satisfied. An example might serve 
to clarify the preceding statements. 


C=0101111111000110 
D=1001001100001010 


is a complementary pair of length 16. Using the interlace method to 
form a pair of length 32 gives 


A=01100011101011111010000001101100 
B=00110110111110101111010100111001 . 


Splitting this A,B pair in the center and adding the quad to this center 


gives 
A*=0110001110101111 11 1010000001101100 
B+=0011011011111010 10 1111010100111001 


Breaking this pair into standard (I I JII IV) form gives 


T=01011111111000110 
II=01010000111001001 


I1I=01011111111000110 
IV=10101111000110110, 


I has 11 ones, Interlacing I with II, whichhas 7 ones, yields 
atotalof 18 ones. Similarly HII has 11 ones, and IV has 10 ones, 
When these are interlaced the code of length 34 has 21 ones, Inthe 
actual computer run the codes of length 16 with 6 ones were converted, 
by complementing, into codes with 10 ones, Therefore boththe A and 
B codes contained 10 ones, An exhaustive list of these and their time 
reversals was generated, Each code of length 16 was then split in the 
center and a zero inserted for the ninth bit, forming an over all code of 
length 1? with 10 ones. A second list was then formed, using a one 
rather thana zero for the inserted ninth bit, forming codes of length 17 
with 11 ones, 

These two lists of numbers were then fed into the program of 


Figure 7,1 at block two, and the regular supplementary check program 
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was used, The codes with 11 ones were used for (I, III) and the 
codes with 10 ones were used for (II, IV). When these codes were 
combined in the interlace scheme, I interlaced with II gave 21 ones 
and III interlaced with IV gave 18 ones. All four interlace combina- 
tions (I II II IV), (IJ I Ill IV), (1 II IV III) and (II I IV III) were 
run exhaustively. The run time was approximately three hours and no 
codes were found, These partial searches for n=34 are documented 
here because they represent two fairly obvious approaches to the prob- 
lem, and it would be very wasteful of time for someone to duplicate this 
effort. The program as given in the appendix can witha slight modifica- 


tion be used for an exhaustive search for kernels of length 34 or fora 


partial search based upon other ‘judicious guesses", 
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CEAr Tt Rees 
CONCLUSIONS 


The study which has been presented here had two principal objec- 
tives. The first objective was to make an exhaustive search for kernels 
of length 26 to determine whether they existed, and if they did exist, 
to determine the number of such kernels, The second objective was to 
develope a methodical scheme of code decomposition from composite 
codes back to their generating kernels, We shall summarize in this 
chapter the ways in which each of the major portions of this dissertation 
is related to the accomplishment of these objectives, Also a number of 
suggestions for further research will be presented. 

The operations group and its family of theorems were originated to 
formalize the elimination of possible redundancies in the search for new 
kernels, The predetermination of the number of ones in each of the 
sequence quadruple vectors was a vital screen that eliminated many 
such redundancies, Also, the starting of both the I and III vectors 
from the same possible supplementary code groups eliminated an ad- 
ditional half of the possible codes, The theorems in the operations 
group chapter also were useful in the code decomposition part of the 
problem, since they show the necessary and sufficient conditions for 
the (I, II) and (III, IV) pairs to be complementary. Theorem 3,4, 
which states that if the sequence quadruple formis (I II I H) then 

(I, II) must be a complementary pair, is most helpful since the process 
of checking for complementarity is very tedious by hand methods, 

Hamming distances were the next major topic considered, They 
were found to be useful in a number of the proofs in the Hamming vector 
chapter, They were alsoa great aid in checking hand decomposition, 
since it was easy to determine whether D(I, IN=DUI1, TV jan / 4 torreon por 


site codes or whether they had the characteristic distance in the case of 


kernels, 111 





Hamming vectors were the most powerful tool developed in this 
paper for code decomposition, The theorems cover all known possible 
decompositions and also give strong support to the conjecture that ker- 
nels of length 50 do not exist. The Hamming vector concept should 
allow a fuller understanding of complementary sequences in general, 
due to its characteristic qualities such as anti-symmetry, and also be- 
cause there are only four possible configurations of Hamming vectors 
for any code pair, Similarly in sequence quadruple form it is easy to 
check whether a decomposed code may be complementary by using its 
Hamming vector. 

Several exhaustive searches for codes of various lengths were con- 
ducted during this research, The most important of these was the 
search for kernels of length 26. This search revealed that only one 
kernel existed, An exhaustive search for codes of length 16 and 20 
revealed that all codes of these lengths were generated from shorter 
codes by standard methods, This gives considerable strength to the 
conjecture that these generating methods are the only ones, 

The possibilities for future research in this field of complementary 
sequences are many. Some of these follow directly from this paper 
while others le quite far afield from the ideas exposed here, An inter- 
esting extension of the present work would be an exhaustive search for 
kernels of length 34. It is felt that the cyclic complementary screen 
which was used in the partial searchis not powerful enough, anda 
better screen is needed before this search can be conducted on an ex- 
haustive basis using today's computers, The writer had considered the 
combination of both the cyclic and supplementary properties as screens 
in the same search, However, a full study was not made of the gain 


which might be obtained by this screening procedure, The programming 
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of such a screen would not prove too difficult and, since the cyclic 
property would logically be the first screen of the two to be used, the 
areas of search could still be catalogued quite easily. 

If it would be possible to determine, before the search, the charac- 
teristic Hamming distances D(I, H]) and D(II], IV) of kernels, these 
would act as a most potent screen in the search for new kernels, It is 
therefore felt that further research on characteristic Hamming dis- 
tances of kernels could prove quite fruitful for future applications, 

It was conjectured that kernels of length 50,130 and 338 donot 
exist, based on several pieces of evidence, which although strong are 
not conclusive, Further investigation along these lines might offer a 
conclusive proof without the necessity for an exhaustive search, since 
an exhaustive search for these length codes approaches the impossible. 
During the investigation of kernels of length 50 the writer noticed an 
oddity which might be worthy of further investigation, 


One kernel of length 10 is 
A=1001010001 
B=10000001 10. 


Divide both the A and B codes into segments of length 2, then 
A=10 01 01 00 01 
B=10 00 00 01 10, 


Now let the symbols 10=A,01=A, 11=B, and 00=B. Forming a code 
of length 50 by this method of using symbols of portions of the quad 
to represent kernels of length 10 gives 

a EON IN 


i ot 


K2=ABBAA. 
This code pair of length 50 satisfies the parity check, and also equa- 
tion 2,5 for the number of ones in both full codes and in the four half 


codes. Several examples were checked in this fashion for both kernels 
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of length 10, and all satisfied the same conditions as did this example, 
but of course none of them satisfied the necessary condition for com 
plementarity. Along with helping to search for codes of length 50, this 
offers a possibility for research in that all kernels seem to have 
various characteristics incommon. However, the quad seems more 
flexible. Of the four known kernels, the quad is the only one thus far 
discovered that can be used to build up composite codes by one method 
of generation and then these codes can be decomposed by a different 
method. Therefore a study of kernel characteristics might lead to 
easier ways of generating new kernels, or at least of determining 
whether they exist for the various possible lengths. 

Another possible field of research is cyclic complementary codes, 
A communications system is more likely to use a cyclic complementary 
carrier modulation instead of straight complementary modulation due to 
the continuous nature of the carrier. Therefore an investigation of 
these codes seems quite in order, 

Two other topics worthy of research in the complementary sequence 
field are the various correlation functions of complementary sequences 
and the frequency spectra of these sequences, For secure communica- 
tions purposes it 1s important to have a code that is noise-like, A 
uniform distribution of ones and zeros gives a noise-like appearance to 
the code. For example, the complementary pair 


A=11010001 
B=11011110 


is a typical complementary code pair of length 8. If the bits of the A 
code are considered ina cyclic fashion one, two, and three bits ata 
time, the results are remarkably uniform as is seenin Table 9.1. 


This uniformity in output would make it extremely difficult for someone 
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monitoring the code to decide if he had noise or an actual signal on his 
receiver, However, it is seenthat the B code is not nearly so uniform 


and would be easier to detect asa signal, 


Peo. 
Bits Number Iwo Bits Number’ Three Bits Number 
] 4 1 2 ba ] 
0 4 10 Z 110 ] 
O01 CZ 101 ] 
00 Z O11 1 
100 l 
010 ] 
001 ] 
000 i 
B code 
] 6 Lt 4 met Z 
0 2 10 Zz 110 Z 
0] Z mop Z 
00 0 O11 Z 
100 0 
010 0 
001 0 
000 0 


Noise uniformity of a complementary 
pair of length 8, 


TAB LE Yee 


An investigation of a large number of codes by the above technique would 
be useful in determining the most noise-like code pairs, where noise- 
like is defined as a uniform output with bits taken one, two, three, etc. 
afatime. lhis might in turn lead to some general characteristmesou 


most-noise-like complementary sequence pairs, = 


ee. 





Another important study would be the autocorrelation function of 
each code pair. This is synonymous with the output of a code detected 
by its matched filter, For example, the autocorrelation of the code A 
given above is 

] 0 -3 0 =! 0 -]1 8 -1 QO -1 OM = 3 1970 1 
where the outputs in time are read from left to right, The autocorrela- 
tion of the B code is 

ai Ooo «60 ] 0 ] 8 1 0 ] O see 0- Ce f 
The autocorrelation functions are important from two standpoints. First, 
if a portion of the system should need to be shut down for repairs or 
routine maintenance, a good autocorrelation function would Alles opera- 
tion on just one carrier, Second, the frequency spectrum of the trans- 
mitted signal is the Fourier transform of this autocorrelation function, 
Therefore, one way to study the spectrum is first to autocorrelate and 
then to take the Fourier transform of the transmitted signals, The more 
widespread the spectrum, the more difficult is a jamming procedure. 

A study of the crosscorrelation functions of complementary sequence 
pairs is also very important because of the difficulty of isolating the A 
code RF carrier from the B code receiver, It was pointed out in 
chapter 4 that because the A code and the B code were orthogonal, at 
exact match cross talk would be no problem. However, when the signals 
are not at exact match cross talk might be a large source of noise, An 
investigation of the crosscorrelation functions of the A and B code 
pairs might lead toa class of complementary sequences where the 
problem of crosstalk is minimized. 

Another class of correlation functions which are important are 
those where cyclic errors are created in the received signal detected by 


the matched filter. These errors could be caused by a linear shift in 
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the phase of the received signal, There may be certain classes of 
complementary sequences where the loss in coding signal-to-noise ratio 
is minimized when noise is caused by an error of this type. This par- 
ticular class of signals would be extremely important in an air search 
radar application, 

A study of the maximum number of different composite codes for 
any length would be an important contribution to the field, because for 
anti- jam reasons it is very important to be able to change codes, and 
the more codes with good characteristics that are available the better, 

The last study to be recommended is an analysis of the autocorrela- 
tion functions, frequency spectra, cross correlation fine none and 
linear phase shift correlation functions for codes of lengths from 100 
to 200, The purpose of this investigation would be to determine if 
certain kernels are capable of forming composite codes of better char- 
acteristics than other’ kernels for radar applications, 

In conclusion, the writer feels that complementary sequence pairs 
will prove very important in future applications to both radar and com- 
munications schemes, However, much research must be done to pick 
optimal classes of these codes for suchapplications, This paper has 
offered many tools to be used by future investigators of complementary 
sequences, These tools have found important applications in their 
present state, but additional investigations along the lines of operations 
groups, Hamming weights and Hamming vectors should further in- 


crease their usefulness, 





GLOSSARY OF SYMBOLS AND TERMS 


(A, B) 
(i: ID 1) 
D(U,V) 


H(U,V) 


AL 
Alter 


Complementary 
code 


Complement a 
code 


Composite code 


Cyclic comple- 
mentary code 


Hamming Distance 
Hamming vector 


Hamming weight 
Kernel 


Quad 


Sequence 
Quadruple 


Supplementary 
code 


Time inverse 


Weight 


A complementary pair of binary sequences, 
see equation 2, 2, 


A complementary pair of binary sequences in 
sequence quadruple form, see equation 3, l. 


The Hamming distance of two binary vectors, 
see equation 4, l, 


The Hamming vector of two binary vectors, 
see the first page of Chapter 5, 


Means interlace two Hamming vectors ina pre- 
scribed manner, see Theorem 5.11, step 2. 


Is transformed into, 


The time inverse of a binary sequence or 
vector A, 


Means to complement every other bit of both 
codes of a complementary code pair. 


A pair of binary sequences which satisfy 
equation 2,2, also called a complementary 
sequence pair, a complementary sequence 
and sometime just a pair or code, 


Means change the sign or take the complement 
of a binary sequence. 


A complementary code whichis reducible to a 
shorter code length by standard methods. 


A pair of binary sequence which satisfy 
equation 6, 3, 


The number of bits in which two binary sequence 
differ, notation is D(A,B), see equation 4.1, 


A binary vector or sequence formed by the 
modulo 2 sum of binary sequences, 


Hamming distance of a binary sequence or vector 
from the null vector, It is the number of ones 
in a binary vector. 


A set of complementary codes which is irreducible, 
The term is sometimes used for a member of the 
set. 


A kernel of length 2, 


A complementary code pair expressed as in 
equation 35.12 


A quadruple of binary sequences which satisfy 
equation 6,1. 


Means to make the first bit of a binary sequence the 
last bit, the second bit next to last, and so on; also 
called time reverse, 
See Hamming weight, 
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EXPLANATION OF FORMAT OF THE APPENDICIES 

Appendix I on matched filters is self-explanatory, Appendix II 
is the result of an exhaustive search for codes of lengh 16. There were 
96 code pairs found by the computer after the program screens had 
eliminated most of the redundancies. The format of Appendix II gives 
the A code, the B code, H(A,B), and H[T, (A,B)] which signifies 
time reversing the A code, It is worthwhile to note that each code pair 
and the pair formed by the time reverse of A were decomposable. 
From these 192 possible decompositions 48 were interlace, 48 were 
time sequence, and 96 were 2" special method, with 48 each of two 
different kinds, 

Appendix III is the result of an exhaustive search for codes of 
length 20. The format of Appendix III gives the A code, the B code, 
if itis required, the time reversal operation used to put H(A,B) in 
standard form, H(A,B) and from which of the 2 kernels of length 10 
the generating codes came. All of the 24 code pairs or their time 
inverse were decomposable. There are 12 from each kernel, 8 of 
each of these interlace and the other 4 are time sequence, 

Appendix IV is a diagrm of all the subgroups of the operations 
group. It starts with the subgroup of order 1 and works upto 
those of order 16. An x inthe row indicates that the operation at 
the top of the column is a member of the subgroup. 

Appendix V is the CDC 1604 computer program in AR format 
and also in machine language for the 26 length search, There are 
slight modifications in the various sections of the program required 
for other length codes, Appendix VI is also in both AR and machine 
language, and is for codes of length 34, This program parallels 


Figure 8.1, whereas the program in Appendix V parallels Fig. 7.1. 
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APPENDIX I 
MATCHED FILTERS 
A Matched Filter is by definition a filter which maximizes the peak 


signal-to-noise ratio. Deterioration of the signal wave form is accepted 


in order to obtain the desired maximum ratio, 

One tool needed for the derivation of the matched filter character- 
istics 1s the Schwarz inequality which is a special case of the Holder 
One representation of the Schwarz inequality is 


Ge 
inequality. 


' 2) etal, 
| x(w) y (w) dw | < | x(w) | dw | y(w) | dw (A 1,1) 
It is to be noted that the equality holds when y(w) is the complex 


conjugate of x(w). Rewriting equation Al.1 gives 


© Zz 
| . x(w) y(w) dw | 


J x(w) [Paw J fy(w)aw 


<a (Ala 


The derivation for the matched filter given here closely parallels 


that of Turin. 28, 


Let f(t) be a signal impressed across a filter whose frequency 
characteristic is H(jw). F(jw), the frequency spectrum of f(t), is 
given by the Fourier transform of f(t). 


F(jw) = } TER ae 


~~ © 


G(jw) is the frequency spectrum of the output after passing through 


the filter 
G(jw) = H(jw) F(jw), = 





or transforming back to the time domain the output signal voltage is 
Jn be 
g(t) = J Hijw) F(jw) eo af. (A 1, 3) 


Since at some time, t= A, g(t) must be a maximum, 


~ 


iw 


max g(t) = g(A) = i” H(jw) F(jw) e 2” af, (A 1, 4) 


To complete the derivation it will be necessary to obtain the total 
noise power, and the total power in the signal, If the noise is assumed 
to be white noise of N watts per cycle, the output noise power density is 

N, = N | H(jw) |°, 


and the output noise power is therefore 


oO 


ly Naf = J Ni [eit bake (A 1.5) 


The total energy contained in the signal is 
& 
= sj Aas dt, 
= © 
which is also 


Geo 5 
E= J |F(jw)|“ae, (A 1. 6) 
by Plancherel's Theorem. It is to be noted that E is a constant not 
dependent upon the filter used. 

The output power signal-to-noise ratio for the maximum output is 


therefore 


ae 
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Dividing this equation by a constant will not effect the time of maximum 


signal output; therefore dividing by E, equation A 1.6, gives 


oO 
re) 


N| g(a) {2 _ [J Fliw) HGw) e IY ae | 2 


a | aw) [far (| F(jwy |? 


{03} = {es) 


df (A 1. 7) 


If the right hand side of equation (A 1.7) is compared.-to the 


Schwarz inequality (A 1.2) with x(w) = F(jw)e we y(w) = H(jw), 


the two expressions are ihe same. But for the expression to be max- 
imum y(w) must equal the complex conjugate of x(w), therefore 
H(jw) = F(-jw) e 7I¥*, 

A matched filter is therefore a filter whose frequency character- 
istic is the complex conjugate of the signal spectrum to whichis matched. 


The output signal for a given signal applied to a matched filter is 


therefore 
OO 


2 (t) 7 | |F (jw) ee jw(t- A) 
“00 lit 
at t= A the output is 
| ee 
g(A) = Jo | F(jw | Gita 


This is the same as equation A 1,6 which is all the energy con- 
taimed in the signal, Theretore at exact match all the energy aneene 


pulse is in the output signal, 
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APPENDIX V 
REMARKS ON COMPUTER PROGRAM 


In Appendix V the addresses from 60000 to 60035 are block 1 
of Figure 7.1, All numbers used in the program are octal. The program 
as written is for codes of length 26, so most of the program deals with 
half length codes of length 13. Words to be changed for different length 
codes are 60000 (mask of ones for half length code as determined by 
weight), 60001 (length of half code minus one [upper half], and length 
of memory word minus length of code [lower half]), 60004 (length of 
half code), 60034 (last sequence in code list). For exactly the same 
reasons respectively change 60015, 60016, 60021 and 60032. 

Addresses 17000 to 17045 separate the sequences into blocks by 
supplementary number. Index register 2 must be entered with 100 
before starting. Addresses to be changed here are 17000 (mask of first 
3 and last 3 bits), 17005 (contains the total number of sequences for 
both half length codes), 17007 (contains the number of sequences with 
7 ones), 17015 (contains a mask for half length codes), 17026 (contains 
the number of sequences with 9 ones), 17032,17033 (have for their 
instruction address the base word 50000 plus the number of half codes 
with 7 ones), 17034 (contains a mask for code length). 

Addresses 10000 to 10106 are the repetitive part of the program. 
The only memory address needing modification here for different length 
codes is 10106 whichcontains SVN a mask equal to the one half code 
length, 

Addresses 7200 to 7252 contain the unpacking subroutine and the 
like-unlike subroutine, 7201, 7207; 72'5, and 7223 all contatnmene 
half code length. Similarly 7204,7207, 7220, 7223 and 7230 all contain 
the full code length or the length modified by one or two to fit the program 


needs, 35 





The 15000 to 15100 series is all possible bit arrangements for 
the first 3 andthe last 3 bits of codes of length 13, while 16000 to 
16100 is a list of supplementary block numbers in one to one correspon- 
dence with this list, 

The 10300 and 10400 series respectively contain the number of 
codes in each supplementary block for the half length codes of weight 
9 and weight 7, 

16000 to 16020 is the final list of the 18 supplementary block 
numbers and is in one to one correspondence withthe 10300 and 10400 


series, 
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